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ABSTRACT 


STRUCTURE  OP  STABLY  STRATIFIED 
TURBULENT  BOUNDARY  LAYER 


The  structure  of  a  stably  stratified  thick  boundary  layer 
developed  in  a  metcro logical  wind  tunnel  is  investigated  experimentally. 
Measurements  of  mean  velocity,  mean  temperature,  turbulent  intensities, 
shear  stress,  heat  fluxes  and  turbulent  spectra  made  at  a  station  78  ft. 
from  the  leading  edge  are  reported.  Turbulent  quantities  were  measured 
by  using  different  hot-wire  probes;  the  measurement  technique  which  is 
a  modification  of  the  procedure  suggested  by  Kovasznay  (1953)  is  des¬ 
cribed.  The  results  show  that  stability  greatly  reduces  the  turbulence 
in  the  boundary  layer. 

The  structure  of  the  wall  layer  is  discussed  in  the  light  of 
Monin  and  Obukhov's  (1954)  similarity  theory  and  Ellison's  (1957)  theory. 
The  results  are  also  compared  with  previous  measurements  in  the  labora¬ 
tory-  and  in  the  surface  layer  of  the  atmosphere  in  stable  conditions. 

It  is  shown  that  mean  flow  and  turbulent  characteristics  of  the  wind 
tunnel  boundary  layer  are  well  described  by  the  similarity  theory,  and 
that  this  theory  provides  a  good  basis  for  the  wind  tunnel  modeling  of 
similar  characteristics  of  the  atmosnheric  surface  layer. 

Measured  spectra  of  lateral  and  vertical  velocity  fluctuations 
arc  shown  to  agree  with  Kolmogorov's  (1941)  similarity  theory  irrespec¬ 
tive  of  the  stability.  The  results  aie  compared  with  Heisenberg's 
(1948)  theory  for  the  equilibrium  spectra.  Spectra  of  temperature 


fluctuations  are  shown  to  have  similar  form  in  agreement  with  Corrsin's 
(1951)  theory.  No  buoyancy  subrange  could  be  identified  in  any  f  the 
spectra. 
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Chapter  I 
INTRODUCTION 

A  fluid  is  considered  to  be  stratified  when  its  mass  density 
varies  with  height.  The  changes  in  fluid  density  may  occur  due  to 
changes  in  temperature,  salinity,  or  some  other  cause.  When  density  in¬ 
creases  with  height,  a  fluid  particle  having  moved  upwards  or  downwards 
from  its  position  of  equilibrium  will  be  subjected  to  a  buoyancy  force 
which  will  try  to  move  it  farther  away.  This  type  of  stratification  is 
called  gravitationally  unstable.  If  density  decreases  with  height,  on 
the  other  hand,  the  buoyancy  force  on  the  shifted  particle  will  try  to 
bring  it  back  to  its  equiliorium  position,  and  the  stratification  in  that 
case  is  called  stable. 

Density  stratified  flows  are  abundantly  met  in  nature,  e.g.,  in 
the  atmosphere,  ocean,  tidal  channels,  ^tc.  When  stratification  is 
caused  by  inhomogeneities  in  the  temperature  field,  it  is  commonly  re¬ 
ferred  to  as  thermal  stratification.  Not  in  all  temperature  varying 
fluid  flows  are  buoyancy  forces  important.  For  example,  in  most  of 
laboratory  heat  transfer  studies,  where  temperature  gradients  are  small, 
stratification  effects  are  found  »o  be  negligible.  A  quantitative 
measure  of  the  effect  of  buoyancy  forces  on  the  flow  regime  is  the 
Richardson  number.  The  Richardson  number  for  the  above-mentioned 
studies  happens  to  be  ve»'v  small.  However,  it  is  not  generally  tlm 
small  for  the  atmospheric  flows  of  meteorological  interest,  in  which 
st  r  it  1  f  1  cat  i  on  effect-  -'list  be  considered.  therefore,  it  j  s  not 


2 


surprising  that  most  of  the  experimental  and  theoretical  work  on 
thermally  stratified  flows  has  been  done  in  connection  with  the  atmos¬ 
pheric  boundary  layer. 

Sufficient  interest  in  measuring  turbulence  in  the  atmosphere 
developed  only  after  1950;  earlier  data  pertain  almost  exclusively  to 
the  mean  velocity  and  temperature  distribution.  Theoretical  develop¬ 
ments  came  even  later.  A  comprehensive  similarity  theory  was  put  forward 
by  Monin  and  Obukhov  (1954),  which  for  the  first  time  provided  a  theo¬ 
retical  framework  for  the  presentation  and  mutual  comparison  of  data  from 
different  sources.  Other  theoretical  models  based  on  the  dynamical 
equations  have  been  proposed  by  Ellison  (1957),  Townsend  (1958),  and 
Monin  (1965),  a.tiong  others.  Sufficiently  accurate  and  simultaneous 
measurements  of  turbulent  intensities  and  fluxes  have  not  been  available, 
however,  to  verify,  conclusively,  the  results  of  these  works  and  various 
assumptions  involved  in  obtaining  them.  It  is  felt  that  laboratory 
studies,  in  which  the  conditions  of  the  model  flows  can  be  artificially 
created,  will,  perhaps,  be  more  useful  than  field  studies  for  checking 
and  supplementing  the  results  of  the  existing  theories,  and  for  develop¬ 
ing  new  ones.  Not  only  can  flow  conditions  in  the  laboratory  be  main¬ 
tained  steadily  for  a  sufficient  time  to  take  accurate  measurements,  but 
also  these  conditions  can  be  systematically  varied  whenever  necessary. 

Only  very  few  laboratory  studies  of  the  stratified  fluid  flows 
exist.  Ellison  and  Turner  (1959,  I960)  investigated  the  entrainment 
characteristics  of  a  mixing  layer  of  salt  solution  in  a  rectangular  pipe, 
hebster  (1964)  studied  the  effect  of  stability  on  The  turbulent  inten¬ 
sities  and  turbulent  fluxes  at  uic  center- line  of  a  wind  tunnel  where 
me  .in  velocity  and  temperature  varied  ltnearlv  His  data,  however, 
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showed  some  peculiarities  of  the  paTticular  wind  tunnel  in  that 
extremely  high  values  of  the  ratio  K../K,.  were  observed,  and  the  flow 
in  the  wind  tunnel  appeared  to  be  highly  developing. 

The  present  study  of  the  thermally  stratified  boundary  layer  in 
the  Fluid  Dynamics  and  Diffusion  Laboratory  at  Colorado  State  University 
has  ben  undertaker  as  part  of  a  long-time  project  which  aims  at  labora¬ 
tory  simulation  of  various  atmospheric  phenomena,  viz.,  mean  velocity  and 
temperature  profiles,  turbulent  intensities  and  scales,  and  turbulent 
diffusion.  The  requirements  of  such  modeling  and  some  previous  experi¬ 
mental  work  in  the  meteorological  wind  tunnel  designed  for  this  purpose 
have  been  reported  by  Cermak,  et  al.,  (1966).  Plate  and  Lin  (1966),  and 
Chuang  and  Cermak  (1966)  shewed  that  the  velocity  profiles  in  the  wind 
tunnel  are  similar  to  those  observed  in  the  atmospheric  surface  layer 
and  are  in  good  agreement  with  Monin  and  Obukhov’s  (1954)  similarity 
theory.  This  earlier  success  in  establishing  a  preliminary  basis  for 
the  modeling  of  mean  velocity  profiles  encouraged  the  present  study  of 
turbulence  structure  in  a  thermally  stratified  boundary  layer,  the  pri¬ 
mary  objective  of  which  is  to  see  whether  a  similar  basis  for  the  model¬ 
ing  of  the  turbulence  structure  also  exists.  It  will  be  shown  that  this 
indeed  is  the  case.  The  results  of  the  present  measurefitepts  will  also 
be  used  to  verify  some  of  the  results  and  assumptions  of  cllison's  (1957) 
theory.  Finally,  the  results  of  the  measured  one-dimensional  spectra 
of  velocity  and  temjier.it  lire  fluctuations  v- i  I  i  be  presented,  .tnd  these 
will  be  compared  with  Kolmogorov '  s  1,1911)  Similarity  theory  for  small 


scale  turbulence. 
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Chapter  II 

THEORETICAL  AND  EXPERIMENTAL  BACKGROUND 

Most  of  the  research  on  thermally  stratified  flows  has  been  done 
in  connection  with  atmospheric  and  oceanographic  studies,  A  brief  re¬ 
view  of  the  various  theoretical  and  experimental  works  is  given  in  this 
chapter.  Some  earlier  measurements  taken  in  the  laboratory  flows  are 
also  discussed.  The  need  of  more  comprehensive  experiments  in  the  labo¬ 
ratory  to  verify  and  improve  upon  the  existing  theories  and  to  develop 
new  ones  is  emphasized.  The  relevancy  of  the  present  experiments  is 
brought  out  in  this  connection. 

2 . 1  Dynamical  Equations: 

The  two  dimensional  boundary  layer  flow  of  a  thermally 
stratified  fluid  is  considered.  The  fluid  is  assumed  to  be  incompress¬ 
ible  although  the  density  is  not  uniform.  This  means  thar  *•:,«;  changes 
in  the  density  are  entirely  due  to  changes  in  temperature,  not  to  changes 
in  the  pressure.  Then,  the  density  and  temperature  fluctuations  are 
simply  related  as 


(2-1) 


in  which  p  is  the  mean  densitv,  T  the  mean  absolute  temperature; 
and  o'  and  t'  are  respectively  the  fluctuations  in  p  >nd  T  . 

Let  U,  V  and  W  denote  the  components  of  the  u.ean  velocity, 
and  u' ,  v*  and  .v '  corresponding  fluctuations  in  x,  y  and  z 
directions.  Then,  the  conventional  boundary  layer  approximation  for  two 
dimensional  me  an  mot  ion  ir.  x,  z- plane  gives  the  following  equations: 
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in  which  v  and  a  are  molecular  diffusivities  of  momentum  and  heat. 

For  turbulent  motion,  following  the  scheme  of  Lumley  and 

Panofsky  (1964),  or  Monin  and  Yaglom  (1965),  equations  can  be  developed 

for  any  arbitrary  number  of  single  point  statistical  moments.  These 

constitute,  however,  an  unclosed  system.  Of  special  importance  are  the 

equations  for  the  second  moments  which  include  the  energy  equation.  With 

an  additional  assumption  of  plar.e-homogeneity  which  is  reasonable  for  a 

well-developed  boundary  layer,  the  following  equations  for  u' ^ ,  v’~^, 

-  ~ 

w'  ,  u'w',  t'  ,  u't'  and  w't'  can  be  obtained. 
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The  dissipation  terms  e  ,  e  ,  etc.,  are  given  by: 
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and  so  on . 

On  the  left-hand  side  of  l.qs .  (2-5)  through  (2-11),  the  first 
term  represents  the  time  rate  of  change  of  the  turbulent  quantity  in 
quest  ion,  which  vanishes  under  steady  state  conditions.  The  various  terms 
in  the  bracket  of  the  second  term  are  due  to  transport  by  diffusion.  Ot 


7 


these,  the  molecular  diffusion  term  is  usually  negligible  except  in 
the  viscous  sublayer  close  to  the  boundary.  On  the  right-hand  side,  the 
first  term  represents  the  effect  of  pressure  fluctuations  in  redistribu¬ 
ting  energy  among  various  fluctuating  components.  This  term  is  absent  in 

______  'Ml  ______ 

Eq.  (2-9).  Terms  like  -u'w'  7—  ,  -w't'  —  ,  etc.,  are  the  so- 

oZ  aZ 

called  production  terms  which  represent  the  work  done  by  the  mean  flow 
against  the  turbulent  shear  stresses.  The  absence  of  production  terms  in 
Eqs.  (2-6)  and  (2-7)  is  noteworthy.  Thus,  the  energy  is  first  supplied 
from  the  mean  flow  to  the  longitudinal  fluctuations  alone;  it  is  then  re¬ 
distributed  to  other  components  through  pressure  fluctuation  terms. 

Terms  containing  the  buoyancy  parameter  ,  are  due  to  thermal 

stratification. 

Equations  (2-5)  through  (2-7),  when  added  together,  give  an  over¬ 
all  turbulent  energy  balance  as 
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in  which  q'“  =  u'“  +  v'*"  +  w'"  ,  and  c  is  the  total  turbulent 

energy  dissipation.  Much  information  on  the  structure  of  turbulence  in 
a  thermally  stratified  boundary  layer  could  be  obtained  if  the  magnitude 
of  various  terms  in  the  dynamical  equations  could  be  determined.  Hut 
such  a  complete  determination  has  not  been  possible  so  far  even  for 
simpler  unstratified  boundary  layer  f lows  Measurement  of  pressure 
fluctuation  terms  has  been  the  greatest  obstacle.  It  is  therefore  more 
convenient  to  sp.nlv  the  total  energy  balance  where  these  terms  add  up  to 
r.ero  because  of  me  continuity  equation  for  the  t  luc  iuat  1  on  motion,  1.0.  , 
1 1 1  'X  ♦  \  '  -v  *  '  •;  such  studies  tor  neutral  boundur.  layer 
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have  been  made  by  Townsend  (1951),  Klebanoff  (1955),  and  Bradshaw  (1967) 
among  others.  It  has  been  suggested  (e.g.,  Rotta  (1962))  that  in  the 
wall  layer  diffusion  terms  are  too  small  and  production  3°^  dissipation 
must  be  approximately  in  balance.  But,  in  the  outer  layer,  diffusion 
terms  might  be  important. 

For  stratified  flows  tiic  additional  term  due  to  buoyancy  is 
present.  The  effect  of  buoyancy  on  the  energy  balance  is  described  by 
the  flux  Richardson  number  defined  as 
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If  one  formally  defines  the  turbulent  diffusivities  of  heat  and 
momentum  as 
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then,  the  relation  between  R_  and  the  ordinary  Richardson  number  R^  , 
is  given  by 
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R-  is  positive  for  stable  stratification  and  negative  for 
unstable  stratification  indicating  that  in  the  former,  energy  is 
drawn  by  the  buovancy  from  the  turbulence  arid  in  the  latter  it  is 
supplied  to  the  turbulence. 


9 


The  energy  budget  in  tin  atmospheric  boundary  layer  from  heights 

of  25  to  100  meters  was  given  by  Panofsky  (1962).  It  was  shown  that  the 
3  ~2 

term  —  (q '  w')  is  import  in t  for  large  negative  Richardson  number, 
but  it  is  quite  small  for  small  Richardson  number.  For  the  near  ground 
layer,  which  is  not  very  strongly  stratified,  Taylor  (1952)  showed  that 
the  production  and  dissipation  are  the  only  important  terms.  Of  course, 
for  strong  stability  or  instability  cases,  the  buoyancy  term  will  be 
important  too.  Then,  Eq.  (2-13)  in  the  absence  of  diffusion  terms  and 
for  steady  conditions  becomes 

c  =  -  ^  ♦  I  k^F  (2-18) 

3  Z  1 

and  can  be  used  to  determine  dissipation  e  ,  which  is  extremely  diffi¬ 
cult  to  measure  directly.  Similarly,  Eq.  (2-9)  can  be  written  as 


e  ,  =  -  2  w't'  (2-19) 

which  can  be  used  to  determine  "dissipation"  et,  ,  of  the  mean  square 
temperature  fluctuations  in  the  region  where  diffusion  terms  are  likely  to 
be  small.  For  a  flat  plate  boundary  layer,  Eqs ,  (2-18)  and  (2-19)  may  be 

f 

expected  to  be  valid  for  j  <_  0,1  ,  excluding,  of  course,  the  viscous 
sublayer.  All  terms  except  c  in  Eqs.  (2-18)  and  ( 2 -  IP 1  have  been  meas¬ 
ured  in  the  present  study.  During  discussion  of  results  the  value  of 
tf ,  as  determined  from  the  measured  spectra  of  temperature  fluctuation 
will  be  compared  to  that  obtained  from  Eq. 
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2 . 2  The  Effect  of  Buoyancy  on  Turbulence: 

Stewart  (1959)  has  examined  the  effect  of  buoyancy  forces  on 
turbulence  by  considering  Eqs .  (2-5)  through  (2-11)  without  the  rela¬ 
tively  unimportant  diffusion  terms.  Under  conditions  of  stable  stratifi- 

~T 

cation,  buoyancy  directly  acts  on  vertical  fluctuations  to  reduce  w' 

as  in  Eq.  (2-7),  which  then  affects,  in  the  same  sense,  -u'w'  and 
2 

u'  respectively,  through  Eqs.  (2-8)  and  (2-5).  Since  the  energy  to 
— 2  — 2  — J 

v'  and  w‘  is  transferred  from  u'  by  pressure  fluctuations,  a 

— I  — 2 

decrease  in  u'  will  be  reflected  in  the  further  reduction  of  w'  and 

2 

also  of  v'  .  The  above  cycle  will  be  repeated  until  a  new  equilibrium 

has  been  reached.  This  simple  mechanism  is  by  no  means  complete  and  has 

to  be  modified  somewhat  by  diffusion  terms.  But,  the  fact  that  the 

effect  of  buoyancy  is  to  suppress  the  turbulence  as  a  whole,  rather  than 

the  vertical  component  aloe:,  is  clearly  brought  out.  Still,  the  effect 

on  the  vertical  component  is  more  significant  than  on  other  components, 

which  causes  anisotropy  under  strongly  stable  conditions. 

Stewart's  (195°'  physical  arguments  about  the  relative  efficiency 

of  the  pressure  fluctuation  mechanism  and  the  decay  mechanism  further 

led  him  to  the  conclusion  that  the  maximum  possible  value  of  under 

stationary  conditions  should  be  considerably  less  than  unity.  This  has 

been  supported  by  all  the  experimental  measurements  of  Rp  although  an 

agreement  on  what  the  maximum  value  should  be,  has  not  been  reached. 

This  maximum  value  of  R^.  has  been  referred  to  as  critical  flux 

Richardson  number  R^-  ,  by  Ellison  (1957)  and  Townsend  (1958), 

c  r 

According  to  Townsend,  the  critical  condition  arises  from  a 
failure  to  achieve  equilibrium  m  the  interactions  between  the 
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temperature  (buoyancy)  field  and  the  turbulent  motion.  He  explains  it 
by  expressing  production,  buoyancy  and  dissipation  terms,  for  given 
gradients  of  mean  temperature  and  velocity,  as  proportional  to  1  ,  ~ 

and  -  powers  of  the  turbulent  intensity,  respectively.  Further,  it 
is  explained  that  introduction  of  radiative  transfer  may  change  this 
critical  condition.  For  a  non-developing  constant  flux  layer  near  a 
restraining  boundary,  Ellison  (1957)  has  estimated  a  value  for  —0.15, 


I  J  cr 

after  making  some  assumptions  regarding  the  ratio  q'  /w'  and  the  rela- 

— 1  — T 

tive  destruction  of  t'  and  q ' ~  in  the  imagined  absence  of  their  pro¬ 
duction  terms.  Experiments  of  Ellison  and  Turner  (1959)  in  the  mixing 
layer  of  salt  and  fresh  water  in  a  pipe  also  have  been  shown  to  favor  a 

value  of  Rf  .15  .  Analysis  of  Proudman  (1953)  of  data  from  the  gulf 
fcr 

of  Kattegat,  on  the  other  hand,  indicated  a  value  of  Rf  =  0.3  . 

rcr 

It  has  been  pointed  out  by  Ellison  (1957)  and  Stewart  (1959)  that 
the  ratio  l^/K^  increases  with  increasing  stability  and  it  becomes 
very  large  near  the  critical  value  of  R^.  so  that,  there  need  not  be 
any  limiting  value  of  R^  .  Reporting  on  some  measurements  of  J.S. 
Turner  in  the  Cavendish  Laboratory,  Townsend  (1958)  has  suggested  that, 
near  critical  conditions,  the  motion  may  take  the  form  of  irregular  gra¬ 
vity  waves  in  which  the  transfer  of  momentum  is  much  more  intense  than 
the  transfer  of  salinity  (or  heat).  Stewart  (1959)  points  out  that  it 
is  doubtful  whether  satisfactory  results  can  be  obtained  by  treating  such 
a  motion  as  turbulence  in  the  ordinary  sense. 
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2.3  Theories  of  Turbulence  in  a  Stratified  Fluid: 


2.3,1  The  Monin-Obukhov  Similarity  Theory: 

Atmospheric  turbulence  data  have  often  been  presented  in 
the  framework  of  the  similarity  theory  of  Monin  and  Obukhov  (1954).  The 
basic  assumptions  are  that  the  flow  is  homogeneous  in  the  flow  direction 
and  fluxes  remain  essentially  constant  with  height.  Then,  consideration 
of  the  pertinent  variables  leads  to  the  following  length,  velocity  and 
temperature  scales : 
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in  which  u*  «  /rQ/p  ,  is  the  so-called  shear  velocity,  and  Hq  is 
the  wall  heat  flux.  The  length  scale  L  is  called  Monin-Obukhov  length. 

According  to  the  Monin-Obukhov  similarity  theory,  all  the  mean 
flow  and  turbulent  quantities  when  nondimensionalized  by  a  proper  combi¬ 
nation  of  L,  u*  and  T„  ,  must  be  universal  functions  of  the  stratifica¬ 
tion  parameter  . 

Thus,  for  the  mean  velocity  and  temperature  profiles,  the  theory 
predicts 
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which  on  integration  give 
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In  Eqs.  (2-23)  and  (2-24),  z^^.  is  a  reference  height  which  can  be 

taken  as  some  appropriate  fraction  of  L  . 

z  z 

Universal  functions  $(?-)  ,  f(— )  ,  etc.,  cannot  be  predicted  by 

L  Lj 

the  similarity  theory  alone.  It  has  not  been  possible  to  determine  the 

form  of  these  functions  from  other  theoretical  considerations.  Only 

2 

their  behavior  in  the  asymptotic  sense  of  ±  °°  has  been  predicted 

(see  Monin  and  Yaglom  (1965)). 

The  so-called  log-linear  law  nas  been  suggested  as  a  first 
approximation  of  the  Eqs.  (2-2?''  and  (2-24)  ,  It  can  be  expressed  as 
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in  which  8  and  8-j.  are  empirical  constants.  There  has  not  been 
much  agreement  in  the  valjes  of  these  constants  as  reported  bv  different 
authors.  A  detailed  account  of  the  results  from  the  atmospheric  data 
has  been  given  by  Monin  and  Yagiom  I19h5).  it  is  noted  that  only  in  a 
very  few  cases  have  i.  ,  u.  ,  and  1,  been  determined  from  direct 
measurement  of  fluxes,  for  stable  conditions  I  ivlor  ;lboO;  and 
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Takeuchi  (1961)  report  values  of  6  ranging  from  2  to  10. 

McVehil's  (1964)  results  indicate  a  value  of  7  .  An  earlier  set  of 
measurements  by  Plate  and  Lin  (1966)  in  the  Army  Wind  Tunnel  at  Colorado 
State  University  also  favor  a  value  of  7  . 

Other  important  results  of  the  similarity  theory  as  applied  to 
the  turbulent  flow  are,  among  others: 


(2-27) 


(2-28) 
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Functions  ,  F?  ,  etc.,  must  be  universal.  Not  many  laboratory 
experiments  have  been  reported  to  check  these  results,  Lllison's  (1962a) 
measurements  in  grid  heated  air  show  onl>  little  effect  of  the  heat  flux 
on  the  ratios  of  the  velocity  fluctuation  components  and  on  the  shear 
stress  correlation  coefficient,  i.ut  strong  effect  on  the  heat  flux  cor¬ 
relation  coefficient.  From  the  atmospheric  boundary  laver.  measure.aent s 
supporting  the  theory  hive  beer,  reported  by  Morin  <  !9(>2)  ,  and  by  ‘Kvun 
and  Yaglom  ;  More  recent  results  of  king  *  196S )  ,  Mordukhovich 
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and  Tsvang  (1966),  and  Cramer  (1967),  who  used  measured  fluxes,  generally 
support  the  similarity  theory.  The  experimental  scatter,  however,  are 
too  large  to  yield  any  well-defined  universal  functions.  In  some  cases, 
height  dependency  is  also  noticeaoie.  It  needs  to  be  shown  that  the 
stability  parameter  ~  also  correlates  turbulence  data  from  laboratory 
flows.  This  is  of  utmost  importance  from  the  point  of  view  of  atmo¬ 
spheric  modeling,  the  requirements  for  which  have  been  discussed  by 
Cermak,  et  al_.  (1966). 

It  will  be  shown  in  Chapter  IV  that  the  present  measurements  of 
turbulence  in  a  stably  stratified  boundary  layer  are  in  good  agreement 
with  the  Mon i n-Ohukhov  similarity  theory.  The  results  will  be  compared 
with  those  obtained  from  the  near  ground  layer  of  the  atmosphere. 

2.5.2  Other  Theoretical  Models: 

Similarity  theories  are  based  on  dimensional  considerations,  and 
as  such  canuot  predict  the  exact  form  of  the  universal  functions.  These 
functions  have  to  be  determined  empirically  from  observed  data,  or  theo¬ 
retically  from  consideration  of  dynamical  equations.  The  latter  approach 
has  been  followed  in  somewhat  different  ways  by  hllison  (1957),  Townsend 
(1958)  and  Monin  (1965).  Diffusion  terms  have  been  neglected  in  all 
these  molds. 

hllison  considered  the  flow  over  an  infinite  rough  plane.  It  is 
assumed  that  the  shear  stress  and  the  heat  flux  in  the  layer  considered 
remain  constant  with  the  height  ;  .  He  considers  the  dynamical 
1  qs .  ( f-‘M  ,  (.'-Hi  and  (1-15*  without  the  dittos;  >n  terns,  and  introduces 

— T  ,  _ 

decay  *!!>•-•  i  j  ,  i.  rid  i  .  tor  t'%  q'-  and  w't'  ,  such,  tn.it  in  the 
absence  of  production  terms  these  quantities  would  begin  to  be  de¬ 
stroyed  at  rates  !  ,  etc.  !n  i?  o,  !,  ,  !  ,  and  S  .  an  defined  a  . 

I  .  s 
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(2-33) 


+  u'w ' 


-  w1 


t'  -  0 


(2-34) 


w't'  ~T  3T  g 

TI  +  w  3z  *  T  * 


(2-35) 


From  Eqs .  (2-33)  through  (2-35)  ,  the  following  expressions  for 
the  ratio  K^/K  and  the  heat  transfer  correlation  coefficients  y  ,  , 
are  obtained. 


q'  w' 


1  -  Rf  (  1  ♦ 

T2  f 

T"  1 1  - 


r  2; 

Rf  1 


3^) 


(2-36) 


Yw '  t  ' 
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-  Rf  1  ♦  J-  4^- 

.2 


2  r~  (  1  -  Rf) 
3 


Further,  assuming  that  I  ,  1.  and  T  remain  in  rough Iv  fixed 

T1 

ratios  even  as  these  vary  with  stability,  and  that  1  and 

*  y 

"" — T  » 

q '  *7  w  1  ‘  5.5  ,  Fllison  obtains  from  Fq  .  ^  -  36 )  a  value  of 


(2-37) 


R.  0.15  .  It  has  been  suggested  that  R might  be  even  smaller 
\r  ■— *  '  hr 


is  the  r  it 


may  increase  near  critical  conditions.  Hut, 


of  i  our  e,  |  i  ,  nv  i  Iso  var\  with  stability  which  can  on  1  v  be  deter  - 
mined  ’  rom  experiments .  Vnothe:  assumption  regard,  nc  the  value  of 
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Tj/T^  in  Eq.  (2-37)  leads  to  the  prediction  of  \tti  "  0.3  ,  in 
the  near-neutral  conditions.  Based  on  the  assertion  that  is  deter¬ 

mined  by  typical  length  and  velocity  scales  of  the  turbulence,  and  tak- 
— 2  h 

ing  (q'  )  for  the  velocity  scale,  Ellison  defines  a  length  scale 

— 7  h 

Lm  =  T2  (q-  )■  .  It  is  predicted  that  1.^  will  be  proportional  to  z 

in  unstable  and  neutral  conditions  and  proportional  to  L  in  very  stable 
conditions.  Another  length  defined  as  L„  =  (t'")^/  ~  ,  will  give 

some  idea  of  the  vertical  distance  traveled  by  fluid  particles  before 
either  returning  towards  their  equilibrium  level,  or  mixing. 

Experimental  verification  of  Ellison's  results  and  the  assumptions 
in  the  model  has  not  been  made,  so  far,  under  proper  conditions  of  the 
model,  Webster  (1964)  discussed  the  theory  in  the  light  of  his  measure¬ 
ments  in  a  stably  stratified  flow  taken  at  the  centerline  of  the  wind 
tunnel.  After  passing  the  air  through  heating  and  shear  grids,  almost 
linear  velocity  and  temperature  profiles  were  obtained  in  the  neighbor¬ 
hood  of  the  centerline.  It  appears  that  Webster's  flow  regime  is  hardly 
the  same  as  visualized  in  Ellison’s  model  which  assumes  a  constant  flux 
layer  near  a  restraining  boundary.  Requirements  of  the  model  flow  have 
been  more  nearly  realized  in  the  present  experiments.  Complete  verifica¬ 
tion  of  the  theory  requires  measurement  of  u'“  ,  v'“  ,  w'*-  ,  t '  “  , 


u'w'  and  w't'  in  the  constant  flux  layer,  in  addition  to  the  mean 
velocity  and  temperature  profiles.  Hus  has  been  accomplished,  and  the 
results  will  be  discussed  in  Chapter  IV. 

Townsend  (1938)  iias  considered  the  turbulent  flow  in  a  stably 
stratified  fluid  far  from  the  restraining  boundaries.  In  contrast  to 
Ellison's  i  193  7)  flow  regime,  Townsend  assumes  his  flow  to  be  homo¬ 
geneous  m  the  direction  oi  shear  and  i nhomogenoous  in  trie  direction 
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of  flow.  The  equations  of  energy  and  mean  square  temperature 
fluctuations  are  considered,  and  dissipation  terms  in  them  are  expressed 
in  terns  of  two  length  scales.  After  introducing  some  simplifying  as¬ 


sumptions,  Townsend  predicts  a  value  of  Rf  s  T  •  Furthermore,  it 

rcr 


1 

7 


is  concluded  that  the  turbulent  intensity  will  be  finite  in  the  critical 
flow,  and  that  a  sudden  collapse  of  turbulent  motion  will  occur  as  the 
limit  is  crossed. 

In  another  theoretical  ir  del.  Monin  (1965a)  assumed  a  flow 
regime  similar  to  that  of  Ellison,  considering  at  the  same  time  the  com¬ 
plete  set  of  dynamical  Eqs.  (2-5)  through  (2-11)  .  A  semi-empirical  hy¬ 
pothesis  is  used  to  express  the  pressure  fluctuation  terms  in  terms  of 
Reynolds  stress,  e.g.. 


P' 


I  fii' 

Uz 


+ 


3w '  j 
3x  / 


u'w' 


(2-38) 


etc.,  in  which  is  a  positive  coefficient.  Monin  obtains  the  fol¬ 

lowing  equations  from  the  solution  of  Eqs.  (2-5)  through  (2-8). 


u't ' 


(2-39) 


It  is  pointed  out  by  Monin  that  u't'  is  probably  too  small,  so  that 
the  right-hand  side  in  Eq.  (2-39)  may  be  approximated  by  two.  However, 
simultaneous  measurements  of  u't'  and  w't'  in  the  atmosphere,  by 
Zubko'ski  and  Tsvang  {196b),  and  hv  Zubkovski  and  Kravchenko  (1967), 
indicate  that  uTfT  is,  in  fact,  several  times  -w't*. 

shows  a  tendency  to  increase  as  stability  increases. 


It  will  be  shown 
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during  discussion  of  the  results  that  the  present  measurements  also 
confirm  the  above  mentioned  results  from  observations.  Therefore,  it 
appears  that  Monin's  assumption  that  u' t'  is  negligible  is  not  jus¬ 
tified  at  least  for  stable  conditions. 

Monin's  use  of  another  hypothesis  that  e,  e  ,  ,  etc.,  are  com- 

— 7  — 7 

pletcly  determined  by  q'  ,  t'  ,  and  the  turbulence  scale  i  ,  leads 
him  to  obtain  an  interesting  relation 

— ,  (2-40) 

B1U*J 

in  which  C  is  a  constant. 


'w't' 


1  +  £■ 


2 ■ ^  Turbulence  Spectra  In  a  Stratified  Fluid: 

In  section  2.1  the  spatial  budget  of  energy  in  a  stratified  shear 
flow  was  discussed.  More  significant  information  on  the  turbulence 
structure  is  obtained  from  consideration  of  the  spectral  energy  budget 
which  describes  the  distribution  of  energy  among  eddies  of  various  sizes. 
While  much  of  the  work  done  in  this  direction  for  unstratified  fluids 
(see  e.g.,  Batchelor  (1953a)  and  Townsend  (1956))  is  relevant  also  for 
stratified  fluids,  added  effects  of  buoyancy  must  be  considered  in  the 
latter  case. 

2.4.1  Spectral  liquations: 

Belgian*'  ll9t>2)  i  i  rst  presented  a  theoretical  model  for  the 
spectrum  of  energy  m  i  plane  homogeneous  flow  of  stably  stratified 
fluid.  In  his  woi>  ,  tiie  etfect  ol  anisotropy  arising  from  the  effect  ot 
buovancy  m  encrg>  removal  iron,  vertical  volocit*  fluctuat  mis  is 
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considered,  and  spectral  equations  of  mechanical  and  thermal  energy  are 
derived.  A  more  detailed  derivation  of  the  spectral  energy  equations 
has  been  given  by  Lumley  (1964b)  who  brings  out  the  complicated  nature 
of  spectral  transfer  terms. 

Following  Phillips  (1965),  the  energy  equations  in  the  wave  num¬ 
ber  space  for  two  dimensional,  steady,  plane  homogeneous  flow  can  be 
written  as 

=  S  (k)  |-^J  -  ~  c  '  (k)  +  f  H(k)  -  2v  k2  E(k)  (2-41) 

^Q-rCk)  3T  3  2 

-  V-H(k)  ¥z  -Tk  £'t  (k)  -  2ak  Vk)  ^ 

in  which  E(k)  is  the  three  dimensional  energy  spectrum  and  ^(k)  the 
corresponding  temperature  fluctuations  spectrum.  In  Eqs .  (2-41)  and 
(2-42)  the  term  on  the  left-hand  side  is  due  to  diffusion  in  the  phys¬ 
ical  space.  On  the  right-hand  side,  the  first  term  is  the  pi-oduction 
term:  S(k)  and  H (k )  being  respectively  the  co-spectrum  of  -u'w'  and 
w't’;  the  secor:d  term  represents  spectral  transfer  of  energy;  and  the 
last  terra  is  viscous  dissipation.  The  appearance  of  an  additional  term 
H(k)  in  Eq,  (2-41),  is  due  to  the  effect  of  buoyancy. 

2.4.2  Form  of  Equilibrium  Spectrum: 

The  so-called  equilibrium  range  of  the  spectrum  is  char¬ 
acterized  by  eddies  of  sizes  much  smaller  than  the  scale  of  energy 
containing  eddies.  According  to  Kolmogorov's  (1941)  theory,  the 
structure  of  such  eddies  remains  largely  unaffected  bv  the  large  scale 


motion.  This  leads  to  a  universal  form  of  the  spectrum  in  the 
equilibrium  range.  Then  Eqs.  (2-41)  and  (2-42)  arc  reduced  to 


-  t'(k)  +  ~r  11  (k)  -  2v  k2  E(k)  =  0  (2-43) 

o  K  i 

-  H(k)  ~  e'T  (k)  -  2uk2  )>T(k)  =  0  •  (2-44) 

The  terms  containing  H(k)have  been  retained  in  the  above  equations 
to  include  any  possible  buoyancy  subrange  in  which  H(k)  is  not  negligi¬ 
ble,  but  S(k)  is  negligible. 

Some  progress  toward  solving  Eqs.  (2-43)  and  (2-44)  can  be  made 
only  after  assuming  some  hypothetical  ^el.  lions  for  spectral  functions 
e'(k)  and  c (k) .  Even  then  it  has  often  been  found  necessary  to 
further  approximate  these  equations  for  certain,  i.uher  vaguely  defined, 
sub-ranges,  and  thus,  determine  the  shape  of  the  spectrum  in  a  piece- 
wise  fashion.  This  requires  the  usual  assumption  of  very  large 
Reynolds  number. 

A.  Inertial  and  Convective  Subranges: 

Most  important  of  the  above-mentioned  spectral  subranges  is  the 
so-ealled  inert i a  1  subrange  in  which  the  energy  is  presumably  trans¬ 
ferred  from  low  to  adjacent  high  ware  numbers  without  anv  loss  or  gain, 
l'he  correspond  ing  ^  T  range  in  the  spectrum  of  temperature  fluctuations 
i s  vailed  "convective  subrange".  According  to  ho  1 mogor  -a, ' s  theory,  a 
t  li  rev  -  d  i  mens  :  ona  1  cnergv  spectrum  m  the  inertial  suhringe  must  have 
t  In  far::' 
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and  the  corresponding  one-dimensional  spectrum  of  u'  -  component,  the 
form 

Fu'  (kl]  =  ai  ,2/3  kl"5/3  (2-46) 

18 

where  A  and  a1  =  A,  must  be  universal  constants.  In  Hq.  (2-46), 

2  TT  f 

is  the  one-dimensional  wave  number  related  to  frequency  as  k^=  -g- , 
The  results  of  the  theory  have  been  verified  by  measurements  of  Grant, 
^tewart,  and  Molliet  (1962)  in  a  tidal  channel;  by  Gibson  (1962)  in  an 
air  jet;  by  Kistler  and  Vrebalovitch  (1961)  in  grid  turbulence;  and  by- 
Pond,  Stewart  and  Burling  (1963)  and  Zubkovski  (1962)  in  the  atmosphere. 
The  value  of  the  constant  a^  ’-  'r  r-een  found  in  the  range  of  0.50  ± 
0.05.  In  fact,  a  -5/3  behavior  has  been  observed  also  in  a  variety  of 
other  measurements  in  the  laboratory  and  in  nature.  The  universality 
of  the  constant  aj^  ,  however,  has  been  questioned  by  Lumley  (1965b)  who 
reported  on  the  measurements  by  Margolis  and  Lumley  (1965)  in  a  curved 
turbulent  mixing  layer  to  show  that  a^  could  deviate  significantly 
from  its  supposedly  universal  value  of  about  0.5  depending  upon  the 
ratio  of  production  to  dissipation.  It  has  also  been  suggested  by  him 
that  similar  deviations  could  occur  under  conditions  of  strong  stability 
and  instability.  In  a  later  article  by  Wyngaard,  Tennekes,  Lumley,  and 
Margolis  (1968),  however,  it  has  been  pointed  out  that  some  serious 
errors  had  occurred  in  measurements  of  Margolis  and  Lumley  (1965),  and 
therefore,  l.umley's  remarks  questioning  the  universality  of  the  constant 
Uj  were  incorrect  . 

Die  present  results  or  one- dimensional  spectra  of  v‘  and  w* 
will  be  presented  m  <h,ipter  !\.  Ihesc  will  be  d  Ci. 'sed  in  the  light 
of  kolmogor.u 's  sunil  nity  the  tv. 
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Corrsin  (1951)  has  extended  Kolmogorov's  ideas  of  local  isotropy 
to  determine  the  form  of  the  temperature  fluctuations  spectrum  in  the 
region  where  convective  and  inertial  subranges  overlap.  It  is  given  by 


*  a  -1/3  t-5/3 

i>T  (k)  =  A^.  et,  e  k 


(2-47) 


in  which  A^,  is  a  universal  constant.  A  value  of  0.35  for  the  constant 
for  the  corresponding  one-dimensional  spectrum  $  ,  (k ^ )  ,  has  been 
given  by  Gibson  and  Schwarz  (1°G5). 

Inertial  and  convective  subranges  will  coincide  only  when  the 
Prandtl  number  is  unity  For  very  small  and  very  large  Prandtl  numbers, 
the  form  of  the  temperature  spectrum  in  the  non-overlapping  range  has 
beer  determined  by  Batchelor  (1959),  and  Batchelor,  Howells,  and 
Townsend  (1959). 


B.  Viscous  Subrange: 

For  wave  numbers  larger  than  those  characterized  by  inertial  or 
convective  subranges,  fqs.  (2-43)  and  (2-44)  reduce  to 


A'  L (k  )  =  0 


(2-48) 


2  . k ~  :  ( k  i  =  (' 


i 2-49) 


Hie  solutions  of  I  qs .  1 7  IS)  and  (2-49)  essentially  depend  on 
the  additional  hypotheses  that  must  be  made  for  st  i  1!  unknown  specti.il 
flux  functions  -  ' ( k )  and  ■  '  iki.  ilassi  •"  ampies  are  those  of 
Obukhoi  (1'iJli,  Heisenberg  il'.MS),  and  Kovaszitav  i  1  *  >  5  S  •  . 
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Heisenberg's  eddy  viscosity  approximation  leads  to  a  behavior  for 

the  spectrum  in  the  large  wave  number  region.  Chandrasekhar  (1949)  has 

generalized  Heisenberg's  treatment  to  obtain  a  form  of  the  spectrum  in 

the  transitional  range,  kovasznay's  (1948)  hypotheses  leads  to  a  form 

of  the  spectrum  with  a  finite  cutoff  wave  number.  A  simple  modification 

of  his  model  by  Pao  (1965)  leads  to  a  spectrum  with  exponential  trail. 

Models  proposed  by  tllison  (1962b)  and  Kraichnan  and  Spiegel  (1962)  also 

have  exponential  high  frequency  ends.  For  the  spectrum,  of  temperature 

fluctuations  Corrsin's  U95i)  extension  of  Heisenberg's  theory  also 
-7 

leads  to  k  behavior,  me  form  of  the  spectrin  in  the  vicinity  of 

I  c5|1/4 

Kolmogorov  wave  number  =  [~|  *  ^oes  not  seem  to  be  very  sen¬ 

sitive  to  the  assumed  relation  for  the  spectral  flux,  and  most  of  the 
proposed  models  represent  the  experimental  data  adequately.  It  will  be 
shown  in  Chapter  IV  that  measured  spectra  of  v',  w',  and  t'  show 
k*7  behavior  as  predicted  by  Heisenberg’s  theory.  For  very  large  wave 
numbers,  however,  validity  of  various  models  is  doubtful. 

C.  buoyancy  Subrange: 

A  buoyancy  subrange  might  precede  the  inertial  subrange  of  the 
spectrum  in  stratified  fluids  under  certain  spc'iil  circumstances. 

For  stable  conditions  two  theories,  one  by  Kolgiano  (1959,  1962)  and 


the  other  due 

to  !  uin  l  ey  (  1 9t>  la  , 

1 9<>5u ) 

,  have  been  proposed. 
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turbulence  or  under  strong  inversions  driven  by  flux  divergence.  In 
steady  shear  flows  as  in  the  present  case,  no  buoyancy  subrange  is 
expected.  Measurements  in  the  atmospheric  surface  layer,  as  well  as 
those  of  the  present  study  to  be  presented  in  Chapter  IV,  confirm  this 
conclusion . 

2.4.3.  Turbulence  Spectra  f rom  the 

Similarity  Theory  of  Monin  and  Obukhov: 

Most  of  the  atmospheric  data  given  by  Russian  workers, 

(e.g.,  Gurvich  (19b0),  Tsvang  (1960),  Zhubkovski  (1962),  etc.,)  on  the 

turbulent  spectra  under  stratified  conditions  have  been  presented  in  the 

light  of  the  Monin-Obukhov  (1954)  similarity  theory.  Accordingly,  the 

dimensionless  spectra  of  various  turbulent  quantities  must  be  some  uni- 

f- 

versal  function  of  dimensionless  frequency  p1  ,  and  stability  parameter 
p  .  This  tneory  should  be  valid  for  the  part  of  the  spectrum  outside 
of  the  dissipation  range  as  against  Kolmogorov's  theory  which  should  be 
valid  for  the  dissipation  range  but  not  for  the  energy  containing  range. 
The  form  of  the  spectrum  in  the  inertial  subrange  is  the  same  for  both. 
Cramer  (I9b7)  reports  that  results  of  similarity  nonr.i  1  i  :.at  ion  for  u' , 
v' ,  and  w'  spectra  in  the  atmosphere  are  quite  successful  for  un¬ 
stable  stratification,  but  unsatisfactory  for  stable  stratification. 

2.5  Ihe  Need  of  laboratory  Stuuies: 

l'he  foregoing  discussion  has  summ.tr  l  led  the  existing  theoretical 
and  similarity  models  for  describing  the  turbulence  in  thermal  1>  stiatt 
tied  t  lows  .  Iliesc  mode  I  s  !»au  1  •.•c»  developed  prim.tr;  !v  for  the  inirpose 
of  study  ing  'he  neat  ground  layer  of  the  atmosphere  ir  ►Inch  tin 
is  sump  t  '  ons  o$  pi  i..--  h  >•••••  ••»,e"e  1 1  v  ot  t  iow  and  constant  fluxes  appear  to 
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bo  fairly  well  realized  if  the  terrain  is  flat  and  uniformly  rough. 
Still,  a  great  numb  r  of  variables  met  under  natural  conditions  make 
it  almost  impossible  to  check  the  finer  points  of  the  theory,  and  to 
provide  some  missing  links,  viz.,  tie  value  of  empirical  constants  or 
the  precise  form  of  the  universal  functions,  etc.  From  this  point  of 
view  alone,  it  is  necessary  to  study  model  flows  in  the  laboratory. 

The  ultimate  aim,  of  course,  would  be  modeling  of  various  atmospheric 
phenomena  once  the  preliminary  bases  for  similarity  between  atmospheric 
and  model  flows  have  been  established.  The  meteorological  wind  tunnel 
at  Colorado  State  University  was  designed  and  buill  for  this  purpose 
in  1962  (see  Plate  and  Cermak  (1963))  and  contitv’^us  efforts  have  been 
made  in  this  direction  since  then  (see  Cermak,  et  al.  ,  (1966)). 

Plate  and  Lin  (1966)  have  shown  that  velocity  distributions  in 
the  wind  tunnel  are  in  good  agreement  with  the  Monin-Obukhov  similarity 
theory;  the  value  of  tne  empirical  constant  6  is  found  to  be  the 
same  as  determined  from  measurements  ir.  the  atmospheric  surface  layer. 
This  oiirprising  success  in  establishing  a  modeling  basis  for  the  mean 
velocity  by  using  1.  3S  length  scale  led  further  to  a  program  (see 
Arya  and  Platt  ,1967))  of  study  of  the  turbulence  structure  in  the 
stratified  boundary  layer.  Hie  presene  study  was  conducted  broadly  ir, 
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Chapter  III 

EXPERIMENTAL  EQUIPMENT  AND  MEASUREMENT  PROCEDURE 

This  chapter  pi  es  a  brief  description  of  the  equipment  used  in 
the  present  study,  and  the  procedure  and  techniques  for  each  set  of 
measurements.  Emphasis  is  given  to  turbulence  ii.casurements .  Sources 
and  nature  of  experimental  errors  are  discussed,  and  corrections 
mentioned,  where  applied.  All  experimental  work  was  done  in  the  Fluid 
Dynamics  and  Diffusion  Laboratory  of  Colorado  State  University. 

3 . 1  Wind  Tunnel : 

Experiments  were  performed  in  the  thick  boundary  layer  of  the 
U.S.  Army  Meteorological  Wind  Tunnel  (Fig.  1)  at  Colorado  State  Univer¬ 
sity.  This  facility  has  been  described  in  detail  by  Plate  and  Cermak 
(1963).  The  boundary  layer  was  developed  along  the  floor  of  the  80  ft 
long  test  section  with  a  6  x  6  ft  cross  section.  At  a  distance  of 
40  ft  from  the  test  section  entrance,  the  test  section  floor  changes 
from  plywood  to  a  40  ft  long  aluminum  plate  which  can  be  heated  or 
cooled  from  he'ow  tt'  .un  temperature  between  2u  and  3;>0°f- .  An  air 
conditioning  system  allows  for  the  ambient  .» i  r  temperature  to  be  main 
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The  adjustable  ceiling  of  the  wind  tunnel  was  set  to  obtain  zero 
pressure  gradient.  The  boundary  layer  was  artificially  tripped  by  a 
saw  tooth  fence  which  was  preceded  by  a  6  ft  section  of  1/2  in.  gravel 
roughness  placed  around  the  perimeter.  This  was  done  in  order  to  ensure 
a  well  developed  boundary  layer  at  the  measurement  station  located  78  ft 
from  the  entrance  of  the  test  section  and  2  it.  before  the  end  of  the 
cooled  plate.  At  this  station,  pitot-static  tube,  thermocouple  and 
hot-wire  p< obes  (Fig  2)  were  mounted  on  a  vertically  movable  carriage, 
allowing  them  to  be  sei  at  er;  neight  above  rhe  floor  by  remote  control. 

3  2  Mean  Velocity  Measurements: 

4  1/8  in.  diameter  standard  pitot-static  tube  was  used  as  a  probe 
for  mca.  velocity  measurements.  It  was  connected  by  flexible  tubing  to 
a  Trans-sonic  Type  1207  Lquibar  pressure  meter  (Fig  3).  This  irstro- 
ment  was  calibrated  against  a  standard  Meriaz.  Model  34  FB  2  TM  micro- 
manometer.  No  detectable  difference  was  observed  between  the  two  within 
the  rar.g<  the  present  >.  -easurements  . 

because  of  its  uisiti  vity  to  dynamic  pressure  changes  due  to 
turbulence,  it  is  very  difficult  to  read  the  Trans- son ics  meter  accu¬ 
rately.  For  better  accuracy,  the  mean  position  of  the  pointer  was 
determined,  for  some  runs,  by  integrating  the  electrical  signal  output 
from  the  meter,  over  .t  period  of  four  minutes,  with  the  help  of  a 
specially  designed  integrator  circuit.  The  circuit  was  first  or  hbr  tted 
bv  putting  the  p  rube  m  the  turbulent  free  sf  teaii 

Non  i  iv  >  os  nii'tj  i  a  p  1 1  oi  -  s  t  at  i  c  tube  measurements  have  been 
discussed  by  tloljst.  >n  {  i \>  38  t  .  Mac  Mi  i  inn  ilboi  and  Satulhorn  I  1  t’tti  •  >  . 

1  veept  '.eta  clioe  to  the  tioor,  most  i, ('  these  errors  arc  negligible . 
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Correction  for  the  turbulence  effects  on  the  pitot-static  tube  measure¬ 
ments  was  found,  following  Goldstein,  about  -2%  at  1/8  in.  from  the 
wall  and  much  less  farther  away,  liinze  (1959)  has  questioned  Goldstein's 
method  which  assumes  lateral-velocity  fluctuations  affecting  the  pitot- 
tube  reading  the  same  way  as  longitudinal  fluctuations,  and  has  sug¬ 
gested  that  the  effect  of  the  former  is  appreciably  less  and  may  be 
even  of  opposite  nature. 

3 . 3  Mean  Temperature  Measurements: 

The  temperature  at  any  point  in  the  thermal  boundary  layer  was 
measured  with  a  copper  constantan  thermocouple  with  its  reference 
junction  in  the  ice-bath.  The  thermocouple  e.m.f.  was  read  out  on  a 
sensitive  mi  1 livoltmeter .  The  temperature  of  the  plate  was  determined 
by  a  set  of  thermocouples  embedded  in  the  plate. 

Errors  in  the  temperature  measurements  were  also  considered. 

Those  due  to  thermal  lag  of  the  thermocouple  were  automatically  elimi¬ 
nated  by  taking  point  by  point  measurements  and  allowing  sufficient 
time  for  the  junction  to  attain  a  steady  temperature.  Following  the 
method  given  by  Kreith  ! 1965} ,  errors  due  to  therms'  radiation  from 

the  junction  to  the  wall  it:  the  presence  of  forced  convection  were 

o 

estimated  and  found  to  he  about  h.S  I-  and  nearly  constant  through  the 
boundary  layer.  No  correction  was  applied  tor  thermal  radiation,  since 
it  would  not  affect  the  relative  temperatures  in  the  boundary  layer. 
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3 . 4  Measurement  of  Turbulence : 

3.4.1  Hot-Wire  Anemometers  and  Related  Instrumentation: 

Four  types  of  hot-wire  probes  were  used  in  the  present 

experiments.  These  were:  a  thin  wire  operated  as  resistance  ther- 

2  2 
morneter  to  measure  1; '  ,  a  normal  wire  to  measure  u't',  and  u'  , 

a  yawed  wire  to  measure  u'w'  and  w't'  ,  and  an  x-wire  for  measuring 

2  T 

v'  and  w'\ 

The  resistance  thermometer  consisted  of  0.000025  in.  diameter, 
0.1  in.  long  wire,  made  of  903)  platinum  and  10%  rhodium,  with  an  elec¬ 
trical  resistance  of  ?bout  1000  ohms  It  was  used  with  a  Wheatstone 
br’dg''  (Fig .  4  j  specially  designed  to  pass  a  small  current  of  0.1  m.a. 
or  less  through  the  wire.  The  high  resistance  of  the  wire  and  the 
small  current  insured  that  the  wire  operated  essentially  cold  and  was 
sensitive  only  to  temperature  fluctuations.  No  detectable  sensitivity 
to  velocity  fluctuations  could  be  found.  A  detailed  description  of  this 
arrangement  has  been  reported  by  Qiao  and  Sandbom  (1964). 

Other  hot-wire  sensors  used  in  the  experiment  consisted  of 
0.0003  in.  diameter,  0.05  in.  long,  80%  platinum  and  20%  rhodium  wires, 
mour^ej  on  ^3  ir. .  diameter  ceramic  probes.  A  two  channel  constant 
temperature  transistorized  anemometer  (Fig.  7)  designed  at  Colorado 
State  University  by  linn  and  Sandbom  (19b7)  was  used.  The  frequency 

response  ot  "  s  ins'  _  t  is  flat  up  to  50,00(1  cps  or  grcatei  . 

lor  measuring  rns  '.allies  of  the  fluctuating  electric  signals, 
a  Bru«' l  and  kie.ir  true  ri.v;  -  vMtmeter  (Fig.  3!  was  used.  The  frequency 
response  of  the  voltmeter  is  t.  it  within  2-200  .*100  cps.  IT  a-  instru¬ 
ment  can  also  be  used  as  i  calibrated  amplifier  w.tii  a  maximum  gain  of 
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Because  of  the  low  level  output  of  the  resistance  thermometer 
wire,  it  was  necessary  to  amplify  the  signal  before  reading  on  the  rms 
meter.  For  this  purpose,  a  Tektronics  Type  122  low-level  Preamplifier 
was  used.  No  amplification  was  needed  for  the  hot-wire  anemometer 
signal  except  when  recording  it  on  magnetic  tape. 

3.4.2  Measurement  Techniques: 

In  a  thermally  stratified  flow,  fluctuations  of  temperature 
are  encountered  in  addition  to  the  fluctuations  of  velocity.  Since  the 
hot-wire  responds  to  both  types  of  fluctuations  at  the  same  time, 
measurement  techniques  are  more  involved  than  with  velocity  fluctuations 
atone.  Methods  of  measuring  turbulence  in  the  presence  of  temperature 
fluctuations  have  been  pointed  out  oy  Corrsin  (1949),  kovasznay  (1953) 
and  Morkovin  (1956). 

In  the  present  study,  the  techniques  developed  earlier  by  Corrsin 
(1949)  and  Kovasznay  (1953)  were  slightly  modified  to  achieve  better 
accuracy.  In  particular,  an  elaborate  direct  calibration  procedure,  to 
be  described  later,  was  adopted  for  determining  wire  sensitivities  to 
velocity  and  temperature  changes,  instead  of  depending  on  theoretical 
relations  which  cannot  properly  take  into  account  all  the  experimental 
conditions.  A  mor>  detailed  description  of  the  measurement  technique 
has  been  given  by  Ary  a  and  Plate  (  19t>,s)  . 

following  the  derivations  given  by  Corrsin  (19491  or  Sandhorn 
(  1967),  one  can  write  the  following  equation  foi  the  iespon.-e  of  a  wire 
held  m  the  \-z  plane  at  an  angle  ••  to  the  direction  of  mean  flow. 
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■)E  ,  1  3E  ,  3E 

—  u  '  +  —  —  w 1  +  —  t ' 

nu  ii  no  5T 


(3-1) 


in  which  e'  represents  the  fluctuation  in  the  wire  output,  and 

3E  9E  ,  9E  .  ,  ....  ,  , 

and  are  respectively  the  sensitivities  of  the  wire  to 

velocity,  yaw  angle,  and  the  difference  in  the  temperatures  of  the  wire 

and  the  air.  In  the  following  discussion  different  wire  configurations 

are  considered  for  the  turbulent  quant ites  that,  can  be  measured. 

A .  Normal  W i r ej 

n  r 

For  a  normal  wire  (  0  =  90°  ),  —  =  0,  it  follows  from  Eq .  (3-1) 

that 


3E  ,  3E 

3U  U  +  3T  1 


or 


e ' 


S  u' 
u 


*  St 


(3-2) 


in  which  symbols  S  and  S  have  been  introduced  for  convenience 
u  t 

After  squaring  Fq.  (3-2)  and  taking  time  averages  one  obtains 


Y  =  u'2  XJ  -  I  u't '  X  +  t 


(3-3) 


in  which 


an  d 


X  - 


t 


(3-4 ) 


I  uuat  ion  (3  3t  i  •»  i  parabola  representing  a  functional  re  l  at  l  en¬ 
voi  t  age  output  of  the  wi . c  ami  the  wire  sensitivity 


snip,  net  ween  t 
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parameter  X.  In  principle,  the  three  unknowns  u'  ,  u't'  and  t' 
can  be  determined  from  Fq .  (3-3)  if  three  pairs  of  X  and  Y  are 
measured.  In  view  of  probable  experimental  errors  involved,  however, 

a  more  than  minimum  data  procedure  suggested  by  Kovasznay  (1953)  was 

2  ~~  ~T 

preferred.  Following  this  procedure  u'  ,  u't'  and  t'  were  determined 

by  the  least-square  method,  from  6  to  8  pairs  of  measured  X  and 

— 2 

Y.  Although  fairly  consistent  results  were  obtained  for  u'  ,  tne  same 

— j- 

was  not  true  for  u't'  and  t'  ,  which  showed  considerable  scatter  with 
often  highly  unlikely  values  in  the  low  turbulence  region.  The  same 
type  of  difficulties  were  earlier  reported  by  Johnson  (1955).  This  led 
to  a  closer  examination  of  the  whole  procedure. 

ihe  fact  that  it  >s  very  difficult  to  make  reliable  measurements 
at  low  overheat  ratios,  which  would  be  necessary  to  obtain  small  X 
values,  offers  an  explanation  for  the  observed  anomaly.  For  X  0, 

the  shape  of  the  parabola  of  Lq .  (3-3)  is  most  sensitive  to  the 

7 

coefficient  of  X  ,  less  sensitive  to  the  coefficient  of  X,  and  least 

"2 

sensitive  to  the  additive  constant  t'“,  especially  when  the  last  two 
coefficients  are  small,  i.e.,  in  the  uder  part  of  the  thermal  boundary 
layer.  Unfortunately,  large  X  values  arc  associated  with  small 

velocities  and  hence  a  low  level  of  turbulence.  Inversely,  if  one  wants 

■  ~s  '  '  f 

to  determine  the  unknown  coefficients  u'~,  u't'  and  t'",  from 
measured  pairs  of  \  end  >  m  the  region  X  0 ,  the  res  hits  would 
he  fairly  reliable  for  u'  ,  less  reliable  for  u't',  and  least  reliable 
tor  t'“.  In  fact ,  it  is  conceivable  that  even  no  it ,  e  values  of  t  ' " 
could  he  obtained  beca.ise  of  small  errors  m  measurements.  Hu-  method 
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u’t'  and  t'"  in  low  speed  and  low  turbulence  level  measurements. 

As  a  proposed  modification  of  this  method,  it  was  considered 

2 

necessary  to  make  an  independent  measurement  of  t'“  by  another  wire 
operated  as  a  resistance  thermometer,  and  to  use  this  information  in 
Eq.  (3-3)  for  determining  the  remaining  two  unknowns  u' ^  and  u't' 
from  normal  wire  measurements.  'If i i s  was  found  to  improve  the  accuracy 
of  results  considerably.  For  this  purpose,  the  two  wires,  one  sensitive 
only  to  temperature  fluctuations  and  the  other  sensitive  to  both  tem¬ 
perature  and  velocity  fluctuations,  were  mounted  parallel  on  a  single 
probe.  The  Eq .  (3-3),  then,  CcJi  be  written  as 


T  ° 

2  =  u'  X"-  -  2  u't'  X 
in  which 


(3-5) 


(3-6) 


2 

which  can  be  solyed  for  u'  and  u't'  by  the  least  square  method. 

For  consistent  results,  4  to  6  independent  measurements  of  X  and 

Y  over  a  reasonably  wide  range  of  overheat  ratios  (1.1  to  1.5)  was 

found  sufficient.  This  method  was  used  for  all  the  measurements  of 
2  "~2 

u'  ,  u’t'  and  t'  reported  in  Chapter  IV. 


B  .  Y awed  Wj 

It  will  be  shown  in  the  following  that  wire  yawed  to 
of  mean  t low  can,  effectively ,  be  used  to  measure  turbulent 
u'w'  and  beat  flux  w't'.  A  yawed  wire  is  xrnsit ivr  tn  n 
t  '  as  in  bq  .  i  s-  1  I  .  'lie  cal  ibrat  ion  procedure  then  would 


the  di 
shear 


u  * 


requi re 


re  c  t i on 
st  res s 

-i  n .  I 


3S 


measurement  of  ^  in  addition  to  ^  and  ^  required  of  a  normal 
wire.  This  is  indeed  very  complicated. 

3F 

Some  simplification  can  be  obtained  if  the  fact  that  depends 

3H 

on  is  considered.  If  the  so-called  "cosine  law"  were  valid,  the 

relation  between  the  two  would  simply  be 


j_  lie 

U  3(7 


3t 

all  coc  0 


(3-7) 


This,  of  course,  is  based  on  the  assumption  that  only  that  component 
of  the  velocity  which  is  normal  to  the  wi-e,  contributes  to  the  heat 
transfer.  That  this  is  not  exactly  so  has  been  shown  by  Sandbom 
(1967),  Webster  (1962),  and  Champagne,  et  al.,  (1967).  It  is  shown  in 
Appendix  A  that  after  allowing  the  parallel  component  of  the  velocity 
to  play  seme  minor  part  in  the  heat  transfer  from  the  wire,  Eq.  (3-7) 
is  modified  only  by  a  constant  factor  as 
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in  which  for  a  given  0  ,  c  is  a  constant  close  to  unity.  An  esti¬ 
mate  for  the  present  experiments  gives  a  value  of  c  =  0.92  for 
0  =  +  45°  (See  Appendix  A) . 

After  substituting  b.q.  (3-8)  in  f.q .  (3-1',  one  obtains 

e'  -  S  (  u ’  ♦  v  w'  cot  o  )  -  St  t'  (3-9) 

o 

from  which  follows,  when  •'  -  4S  , 

e '  -  S  (  u*  ♦  c  w*  ) 
u 


S  t 


(  3-  KM 
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and,  when  -  -  45°, 

e'  -  S  (  u'  -  c  w'  )  -  S  t '  (3-il) 

2  u  t 

Squaring  and  time-averaging  Eqs .  (3-10)  and  (3-11)  and  subtract¬ 
ing  the  former  from  the  latter,  will  result  in  the  following  equation: 

Z,  =  -  4  c  u'w*  X“  +  4  c  w't'  X  (3-12) 


in  which 


and 


X  = 


(3-13) 


Liouation  (3-12)  is  similar  to  Hq.  (3-5)  and  can  be  solved  to 
determine  u'w'  and  w't'  by  the  least  square  method  from,  say,  6  to 
10  measured  pairs  of  X  and  V,  All  calculations  were  done  with  a 
IBM-6400  digital  computer.  The  computer  program  for  L:q .  (3-5)  or 
(3-12)  is  given  in  Appendix  B. 

C.  X  -  Wires: 

I'he  vertical  ami  lateral  components  of  the  velocity  fluctuation 
are  measured  by  ope  rat  mg  two  we  1 1 -mat cned  wires  mounted  on  a  single 
probe  in  the  shape  of  X  as  ,n  conventional  hot-wire  anemometry. 
Assuming  t  bat  *  he  two  k  :  res.  have  the  >vun  i  i  i  wi  les.  it  toi  lows 

f  rom  i  qs  .  ;  "i  ■»  ■  an  a  i  •>  1 1 1 1  1 h a { 
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e'  -  e'  -  2  c  S  w' 


(3-14) 


or 


\  w'“ 


2  c  S 


(3-13) 


in  which  is  the  rms  of  the  instantaneous  difference  of  the  sig¬ 

nals  from  the  two  wires, 

.r=f 

Similarly,  Vv'  is  measured  by  operating  the  X-wires  in  the 
x,y-  plane.  The  method  is  simple,  but  it  is  necessary  to  insure  that 
the  sensitivities  of  the  two  wires  do  net  differ  by.  say,  more  than  5"° 


3.4.3  Calibration  Procedure: 

The  calibration  procedure  for  the  hot-wire  in  the  present 

8F  8E 

experiments  was  designed  to  furnish  values  of  and  --  in  the 
range  of  velocities  and  temperatures  which  the  wire  encountered  during 
actual  measurements.  For  this  purpose,  a  set  of  calibration  curves 
was  obtained  by  placing  the  wire  in  the  turbulent  free  stream  outside 
the  boundary  layer  and  measuring  b  vs  U  for  various  wire  tempera¬ 
tures,  keeping  the  air  at  a  constant  temperature.  From  these  curves, 


•IF 

Jli 


was  determined  by  measuring  slopes  off  th-'  F  vs  1!  curves,  or 


by  differentiating  a  mathematical  curve  based  on  king’s  law  nr  a 
similar  power  law  relation,  fitted  through  the  calibration  points.  Hie 
latter  approach  has  been  found  to  be  more  convenient  and  to  give  better 
and  more  consistent  results. 

Ine  sensitivity  ^  ,  was  similarly  measured  after  rep  lot  ting  the 

elibration  d  it  a  in  the  form  of  i  vs  (T  -  I  )  for  constant  U  values 

w  a 

I'he  results  of  a  typic.il  c  a  1  i  b  i  at  i  e.i,  of  a  O.OOOt  inch  diameter  pi  it  i  nun.- 
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rhodium  w i  re  arc  shown  in  F;igs.  6,  7,  and  8.  The  calibration  procedure 
was  repeated  for  each  new  wire  and  also  when  the  wire  calibration  was 
found  changed  during  roucine  checks,  which  were  performed  every  day. 

3.4.4  Errors  in  Turbulence  Measurements : 

It  is  known  that  hot-wire  measurements  of  turbulence  are 
subject  to  various  types  o:r  errors.  Those  of  random  nature  occur  in 
the  process  of  calibration,  reading  instruments,  alignment  of  probes, 
etc.  Systematic  errors  are  introduced  due  to  finite  wire  length, 
gradients  of  velocity  and  temperature,  proximity  of  solid  walls,  high 
intensity  of  turbulence,  inadequate  frequency  response  of  instruments, 
and  due  to  approximations  involved  in  evaluation  and  calibration  tech¬ 
niques.  Most  of  these  errors  become  large  only  when  measurements  are 
made  very  close  to  the  floor,  say,  inside  the  viscous  sublayer.  The 
present  study  pertains  to  the  region  outside  this  thin  layer.  Errors 
due  to  gradients  of  velocity,  temperature,  and  turbulent  intensity  are 
important  only  in  measurements  with  yawed  and  X-wires.  X-wire 
measurements  are  furthermore  affected  by  the  correlation  error  due  to 
finite  separation  distance  between  twe  wires.  For  these  reasons,  yawed 
wire  and  X-wire  results  are  subjected  to  a  greater  uncertainty  level 
than  those  of  normal  wires.  In  the  absence  of  any  suitable  technique 
avail  tble  for  correcting  not -wire  measurements  in  combined  temperature- 
velocity  flow  fields,  no  corrections  were  applied  in  the  present  tur¬ 


bulence  measurements. 
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3.5  Measurement  of  Frequency  Spectra: 

Frequency  spectra  of  lateral  and  vertical  velocity  fluctuations 
and  temperature  fluctuations  were  measured  at  selected  points  in  the 
boundary  layer.  For  this  purpose  fluctuating  signals  were  recorded  on 
FM  magnetic  tape  after  proper  amplification,  and  were  later  analyzed  on 
a  Bruel  and  Kjear  Type  2109  spectrum  analyzer  (Fig.  3).  This  instrument 
const  ts  essentially  of  a  set  of  passive  filters  in  the  range  of  16  to 
10,000  cps.  The  filters  are  of  octave  type  varying  in  band  width 
approximately  proportional  to  the  central  frequency. 

The  spectral  measurements  are  subjected  to  errors  due  to  finite 
wire  length,  filter  band  width,  and  noise.  Due  to  finite  length  of  the 
wire,  two  or  more  eddies  of  small  size  may  strike  the  wire  simultuneously 
with  the  result  that  the  measured  power  would  be  larger  than  the  true 
power.  For  this  error  to  be  small,  wire  length  should  be  small  compared 
to  the  wave  length  of  a  particular  wave-  number.  The  present  measure¬ 
ments  do  not  cover  wave  numbers  larger  than  1100  ft  1  which  is  well 

» 

below  the  value  of  1500  it  corresponding  to  a  wave  length  equal  to 
the  wire  length  (0.05  inch).  No  corrections  have  been  applied  o  the 
present  spectral  measurements. 


RESULTS  AND  DISCUSSION 


The  results  of  the  measurements  in  a  stably  stratified  boundary 
layer  are  presented  in  this  chapter.  Some  broad  characteristics  of 
the  momentum  and  thermal  boundary  layers  in  the  Colorado  State 
University  Army  Meteorological  Wind  Tunnel  which  was  used  for  this  study 
are  discussed.  The  turbulent  intensities  in  the  stratified  boundary 
layer  are  compared  with  csose  obtained  under  neutral  conditions. 
Measurements  of  turbulent  fluxes  and  correlation  coefficients  in  the 
lower  half  of  the  boundary  layer  are  presented,  and  the  effect  of 
stability  on  these  quantities  is  discussed.  Probable  errors  in  the 
measurements  are  also  mentioned. 

Emphasis  is  given  to  the  so-called  wall  region  of  the  boundary 
layer.  The  results  of  the  mean  velocity  and  temperature  profiles, 

KH 

turbulent  intensities  and  fluxes,  and  other  parameters  such  as  -r—  , 

M 

R  ,  R^,  etc.,  are  presented  as  functions  of  the  stability  parameter 
.  These  are  coxpared  with  the  Monin-Obukhov  similarity  theory  and 
also  with  other  measurements  in  the  laboratory  as  well  as  in  the 
atmospheric  surface  layer.  The  nature  of  various  universal  functions 
is  brought  out.  The  theoretical  model  by  Ellison  (1957)  is  also  dis¬ 
cussed  in  the  light  of  present  measurements.  Finally,  measured  one- 
diraensior.al  spectra  of  v*,  w',  and  t'  are  reported;  these  are 
discussed  in  the  light  of  Kolmogorov's  similarity  theory. 

Various  parameters  of  the  boundary  layer  have  been  summarized 


in  Table  I . 
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4 . 1  Broad  Characterits ics  of  the  Boundary  Layer: 

A  striking  feature  of  the  boundary  layer  which  develops  in  the 
Army  Meteorological  Wind  Tunnel  is  its  large  thickness  (more  than  24 
in.  at  the  measuring  station).  This  makes  it  specially  suitable  for 
obtaining  large  Richardson  numbers,  both  negative  and  positive,  and 
hence  for  studying  the  effect  of  thermal  (density)  stratification  on 
the  turbulence  structure. 

In  most  theoretical  studies,  the  assumption  of  plane  homogeneity 
of  the  flow  is  made.  In  the  true  sense  such  a  condition  cannot  be 
achieved  in  a  flat  plate  boundary  layer  which,  because  of  the  shearing 
stress  at  the  wall,  must  be  increasing  with  distance.  But  at  a 
sufficiently  large  distance  from  the  leading  edge,  velocity  and  tem¬ 
perature  profiles  and  also  turbulent  intensities  no  longer  change 
noticeably  with  distance.  Thus,  for  all  practical  purposes,  plane 
homogeneity  can  be  assumed  to  exist  at  long  distances.  For  this  reason, 
the  present  measurements  were  made  at  a  station  where  the  momentum 
boundary  layer  has  developed  for  78  ft,  a  J  the  thermal  boundary  layer 
for  38  ft  length  as  shown  in  Fig.  9.  The  velocity  and  temperature 
profiles  at  this  station  are  compared  in  Figs.  (10)  and  (11)  with  those 
at  a  station  8  ft  further  upstream  to  show  that  the  difference  is  .ndeed 
very  small  and  that  plane  homogeneity  can  be  assumed  to  have  been 
realized  in  the  lower  6  in.  of  the  boundary  layer. 

Two  dimensionality  of  the  flow  near  the  center- line  of  the  wind 
tunnel  was  checked  by  taking  horizontal  traverses  across  the  wind 
tunnel  at  several  heights.  Deviations  from  the  center-line  velocity  of 
not  more  than  1%  were  observed  within  6  in.  on  both  sides  of  the  center- 
line.  There  arc,  however,  indications  of  large  comer  effects  and 
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secondary  circulation  at  low  wind  speeds.  Therefore,  ambien'  velocities 
smaller  than  10  fps  were  not  used  even  though  such  velocities  would  have 
been  desirable  for  obtaining  a  still  larger  range  of  stabilities  than 
covered  by  the  present  experiments. 


4.1.1  Distribution  of  Ri chardson  Numbers: 

A  quantitative  measure  of  the  thermal  stability  is  the 
Richardson  number.  As  defined  in  Eq.  (2-17),  it  is  a  local  parameter. 
For  a  given  temperature  difference  between  the  wall  and  the  ambient  air, 
a  desired  range  of  can  be  obtaine<  by  a  proper  choice  of  ambient 

air  velocities.  In  the  present  :tudy  velocities  of  approximately  30, 
20,  and  10  fps  were  used  to  cover  a  range  of  stabilities  from  near 
neutral  to  moderately  stable.  The  distribution  of  the  in  the 

boundary  layer  for  these  three  cases  is  shown  in  Fig.  12.  A  similar 
plot  on  a  linear  scale  showed  that  for  small  z,  R^  varies  linearly 
with  z  just  as  has  been  observed  by  Plate  and  Lin  (1966),  and  from 
measurements  in  the  atmosphere  as  reported  by  Lumley  and  Panofsky 
(1964).  Batchelor  (1953b)  showed  that  if  suitably  defined,  the 
Richardson  number  can  be  a  reference  parameter  for  the  whole  flow 
field  and  not  merely  a  local  quantity.  Following  Ellison  and  Turner 
(1960),  one  can  define  an  overall  Richardson  number  for  a  layer  of  say 
o  <  z  <  h  as 


R. 
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g. 

T 


T  )h 
o 


(4-1) 


in  which  T  is  the  average  absolute  temperature  in  the  layer.  Thus, 

3 

for  the  whole  boundary  layer. 
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(4-21 


can  be  used  as  .*»  reference  parameter.  Values  of  R.  for  three  cases 

represented  ir-  Fig.  12  are  indicated  tn  the  same  figure  for  comparison 

with  the  local  R.  values.  R.  or  R.  can  be  a  useful  parameter 
1  xh  16 

in  the  outer  region  of  the  boundary  layer  where  flow  characterstics  do 
not  depend  only  on  a  local  parameter  such  as  R^  . 

4.1.2  Mean  Velocity  and  Temperature  Profiles: 

Tne  mean  velocity  and  temperature  profiles  in  the  bound¬ 
ary  layer  for  three  ambient  velocities  are  shown  in  Figs.  (13)  and  (14). 
In  an  earlier  study  using  the  same  wind  tunnel,  Tillman  (1967)  observed 
that  in  the  neutral  case  velocity  profiles  for  different  ambient 

velocities  (20  to  40  fps  in  his  case)  become  similar  when  plotted  as 
U  z 

—  vs  -j-  .  However,  this  type  of  similarity  is  not  achieved  by  the 

CO 

velocity  profiles  in  stable  conditions  as  shown  in  Fig  .  13.  Especially 
the  temperature  profiles  show  some  peculiar  behavior  near  the  edge  of 
the  thermal  layer.  Here,  the  temperature  gradient  instead  of  decreasing 
monotonically  to  approach  zero,  shows  some  sudden  increase  before  finally 
leveling  off.  This  is  perhaps  due  to  imcomplete  mixing  of  the  air  in 
the  core  region  of  the  wind  tunnel,  when  a  residual  stratification  may 
exist  even  after  recirculation.  For  this  reason,  thickness  of  the 


thermal  layer  6,^.  could  not  be  determined  precisely.  It  was  determined 
approximately  after  correcting  for  the  observed  defect  in  the  tempera¬ 
ture  profile  near  the  »'d  e  of  the  thermal  layer.  For  this,  the  normal 
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boundary  layer  type  of  temperature  profile  which  was  observed  over  most 
of  the  thermal  layer  was  slightly  modified  in  the  outermost  region  so 
that  the  temperature  gradient  decreased  to  zero  monotonically .  In  this 
way,  a  value  of  6^.  =■  0.65  6  was  obtained  which  may  be  in  error  up  to 
10%.  Mean  velocity  and  temperature  data  from  the  wind  tunnel  together 
with  other  flow  parameters  for  different  runs  are  tabulated  in  Table  II. 

4.2  Turbulence  Characterstics  of  the  Boundary  Layer: 

4.2.1  Turbulent  Intensities 

The  r.m.s.  velocity  and  temperature  fluctuations  were 
measured  following  the  procedure  described  in  Chapter  III.  These  data 
are  given  in  Table  III.  The  distribution  in  the  boundary  layer  of 
these  fluctuations  normalized  with  respect  to  U  and  (T  -  T  )  are 

*  oo  oo 

shown  in  Figs.  15,  16,  and  17.  Also  represented  in  the  same  igures 
are  for  comparison,  some  results  of  measurements  obtained  under  neutral 
conditions  which  were  made  with  the  same  probes  and  instrumentation  set 
up.  The  following  discussion  is  made  without  any  particular  reference 
to  the  wall  layer  and  various  theories  proposed  for  that  region.  The 
structure  of  the  wall  layer  will  be  taken  up  separately. 

Even  in  the  absence  of  density  stratification,  turbulent  inten¬ 
sities  are  known  to  vary  not  only  with  the  height  z  from  the  wall,  but 

Uoo  6 

also  with  the  characteristic  Reynolds  number  R  =  -  ,  of  the 

e,  v 

0 

boundary  layer.  The  latter  effect,  though  insignificant  when  R  is 

e6 

large,  becomes  important  for  small  R  .  Under  conditions  of  stable 

5  5 

stratification,  it  is  seen  from  Figs.  15,  16,  and  17  that,  as  l)  is 

ao 

decreased,  the  effect  of  both  increased  and  decreased  R  ,  as 

*6  C6 


45 


a  consequence,  is  to  reduce  the  turbulent  intensities.  The  effect  of 

stability  alone  can  be  seen  when  turbulent  intensities  for  the  stable 

and  neutral  cases  are  compared  for  the  same  Rg  in  two  cases.  It 

6 

can  be  seen  that  for  =  30  fps  there  is  hardly  any  change,  but  it 
becomes  more  and  more  significant  as  the  velocity  is  reduced  to  20  fps 
and  then  to  10  fps.  Both  lateral  and  vertical  components  seem  to  be 
equally  affected,  in  spite  of  the  fact  that  the  bouyancy  acts  directly 
only  on  the  vertical  fluctuations.  The  reason  for  this  has  been  given 
in  section  2.2.  The  change  in 
than  the  change  ir 


v^/U 


\T- 


•  /U  with  z/6  is  more  significant 


and 


\R‘ 


•  /U  which  for  R.  =  0.025 


and  0.108  are  seen  to  be  almost  constant  in  the  region  0.05  <  z/6  <  0.5. 

The  temperature  flucuations,  in  contrast,  decrease  rapidly  with  z/6 

and  do  not  vary  significantly  with  R.  .  Perhaps,  the  explanation  for 

*6 

this  lack  of  variation  with  stability  is  that  the  temperature  fluctu¬ 
ations  are  largely  dependent  on  the  temperature  gradients,  which  are  not 
significantly  different  in  the  three  cases,  where  the  difference  in 
temperature  of  the  plate  and  the  ambient  air  is  almost  the  same. 

The  longitudinal  fluctuation  intensities  as  normalized  with  the 
local  velocities  are  shown  in  Fig.  18.  A  sharp  increase  in 
towards  the  wall  is  noticed  for  all  three  stability  cases  and  at  a 


u'  /U 


given  distance  it  is  seen  to  decrease  as  R.  increases.  The  dis* ri- 

1 6 

bution  of  the  turbulent  energy  in  the  boundary  layer  is  shown  in  Fig. 
19,  which  clearly  demonstrates  the  effect  of  stability  on  the  turbulent 


energy. 
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which  was  taken  to  be  the  value  of  the  wall  flux,  "he  thickness  of  this 

layer  is  found  to  vary  from  about  0.5  in  to  2  in.  with  a  tendency  to 

increase  with  stability.  Both  u’w '  c;id  w't'  are  seen  to  be  affected 

considerably  by  the  stability,  fart,  of  the  'hange,  of  course,  might  be 

due  to  change  in  R  . 

e6 

The  heat  flux  in  the  direction  of  flow  is,  representeu  u<  Fig.  24. 

It  seems  to  be  continuously  increasing  as  the  wall  is  approached  where 

the  magnitude  is  several  times  that  of  vertical  flux  as  shown  in  Fig. 

25.  The  quantity  u't'  drops  sharply  with  distance  from  the  wall  and 

it  approaches  a  zero  value  towards  the  edge  of  the  thermal  boundary 

layer.  The  effect  of  stability  on  u't'  is  much  less  than  that  on 

w't'  .  Measurements  of  Webster  (1964)  in  the  wind  tunnel  and  those  of 

Zubkovski  and  Tsvang  (1966)  and  Zubkovski  and  Kravchenko  (1967)  in  the 

atmosphere  also  show  that  u't'  is  several  times  -w't'  under  stable 

u  *  t  * 

conditions  and  that  the  ratio  -  ,  ■■  ,  increases  with  an  increase  in  the 

w't 

stability.  A  physical  explanation  of  this  has  not  been  given.  But, 
comparison  of  the  production  terms  in  f.qs.  (?-!0i  and  (2-11)  indicates 
that  j't1  is  produced  from  che  work  done  by  the  mean  flow  against 

be.;  -u'w'  and  -w't'  ,  while  w't'  is  produced  only  from  the  work 

2  _ 

done  against  w!  .  The  sign  of  the  production  terms  for  u't’  and 

w't*  also  explains  why  the  two  fluxes  should  be  of  opposite  signs. 

The  correlation  coefficients 
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Vf 


,  as  calculated 


from  the  flux  and  intensity  measurements  aie  si. own  in  Figs.  26,  27.  and 
28.  In  general  -y  ,w,  and  ”Yw't'  s*low  a  tendency  to  decrease  with 
the  increase  in  the  stability.  Their  variation  with  z/6  also  becomes 
more  conspicuous  as  the  flow  becomes  more  stable. 

Some  phenomenological  theories  use  eddy  diffusivities  and 

as  defined  in  Eqs.  (2-15)  and  (2-16).  It  has  often  been  found 
necessary  to  make  an  assumption  for  the  magnitude  of  the  ratio 
and  its  variability  with  height  and  stability.  Eddy  viscosity  in  the 

Si  s 

normalized  form,  -Q-y  .  and  the  ratio  jr—  are  presented  in  Figs.  29 

CD 

and  30.  In  Fig.  29.  near  neutral  case  of  U  *  30  fps  (R  =  3.02 

00  r  ex 

5  6 

x  10  ) ,  is  to  be  compared  with  the  results  of  Bradshaw  (1967)  for  the 

neutral  boundary  layer.  The  agreement  is  very  close.  The  effect  of 

stability  in  reducing  the  turbulent  transport  of  momentum  is  clear’ y 

S 

seen  from  Fig.  29.  Fig.  30  shows  that  the  ratio  —  decreases  with 

Si  ^ 

stability.  The  variation  of  -jr-  with  the  distance  from  the  wall  is 

s 

seen  to  become  large  only  near  the  wall.  In  most  of  the  inner  half 
region  of  the  boundary  layer,  its  magnitude  is  about  constant  and  equal 
to  0.85  for  R.  *  0.0112,  0.70  for  R.  »  0.0247,  and  0.60  for 

s, 

R.  »  0.0883.  The  effect  of  thermal  stratification  on  the  ratio  - — 

S  Si 

will  be  discussed  further  in  connection  with  the  structure  of  the  wall 
layer . 
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4.2.4  Probable  Errors  in  the  Flux  Measurements: 

Different  types  of  errors  that  are  common  in  hot  wire 
measurements  have  been  mentioned  in  Chapter  ill.  In  view  of  the  fact 
that  the  results  wore  obtained  after  a  complicated  measurement,  cali¬ 
bration  and  data  reduction  procedure,  it  is  difficult  to  make  an  ac¬ 
curate  analysis  of  the  errors  in  the  flux  measurements.  It  is  expected 
that  th:  use  of  more  than  the  minimum  number  of  data  that  was  necessary 
for  obtaining  the  final  results  reduces  the  random  errors.  Moreover, 
the  detailed  calibration  of  the  wire  at  different  temperatures  provides 
consistency  checks  on  the  measured  sensitivities  of  the  wire.  Some 
independent  checks  on  the  ac  u racy  of  the  shear  stress  measurements 
were  made  by  comparing  values  of  u„  as  obtained  from  hot-wire  measure¬ 
ments  with  the  values  of  u*  inferred  from  the  "wall  layer"  velocity 
profiles.  This  was  possible  only  fo  »  30  and  20  fps  in  which  case 
a  log-law  could  be  assumed  to  be  valid  for  small  z  .  The  measured 
values  of  u„  weie  found  to  be  in  good  agreement  with  those  inferred 
from  velocity  profiles  as  seen  in  the  following  Table. 

Ti.le  (4-1) 

U  u,  from  m.asure-  u„  from  Mean  Vek.  it,'  n  calculated 


fps 

ments  s  u’w’, 

•*ps 

Profiles, 

fps 

from  Ludwieg  and 
Tillmann’s  For,  la, 
fps 

30 . 0 

o 

o 

0.97 

0.967 

20.0 

0.584 

0.57 

0.656 

10. 7 

0 . 240 

_ 

0.337 

su 


Furthermore,  for  the  ne*r  neutral  case  of  =  50  fps,  u*  in  the 
present  case  compares  very  well  with  the  value  of  u*  calculated  from 
Ludwieg  and  Tillmann's  (19S0)  formula. 

These  checks  indicate  that  the  probable  errors  in  the  present 
measurements  of  u'w'  may  not  be  significantly  greater  than  those 
common  in  the  corresponding  measurements  in  the  neutral  flow.  Errors 
in  w't'  would  be  somewhat  greater  than  those  in  u'w*.  As  an  inde¬ 
pendent  check  on  the  accuracy  of  heat  flux  measurements,  the  wall  heat 
flux  can  also  be  calculated  by  usxng  the  following  equation  given  by 
Reynolds,  Kays  and  Kline  (1958)  for  heat  transfer  from  boundaries  with 
step  discontinuity  of  temperature 


in  which  —  is  the  exponent  in  the  power- law  approximation  of  the 
m 

velocity  and  temperature  distributions  in  the  boundary  layer:  a  value 
of  it  =  6  was  observed  in  the  wind  tunnel  by  Plate  and  Lin  (1966). 
Equation  (4-3)  is  based  on  the  assimptions  that  velocity  and  temperature 

Si  *h 

profiles  are  similar,  and  —  =  1  •  Since  the  ratio  in  the  present 

Si  Si 

case  is  consistently  smaller  than  i.nity,  its  actual  average  value  was 
considered,  and  Eq.  (4-5)  modified  as 
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was  used  for  calculation.  In  view  of  the  observed  peculiarity  of 

temperature  profiles  near  the  edge  of  the  thermal  layer,  the  thickness 

‘■'T 

of  the  thermal  layer  6^  ,  and  consequently  the  ratio  ,  could  not 
be  determined  accurately.  But  Eq.  (4-4)  is  hardly  sensitive  to  ^T/6  > 
sir.ce  the  later  is  raised  to  a  power  of  about  -  i  .  A  /alue  of 

5 

T/6  -  0.65  as  discussed  in  section  4.1  has  been  used.  An  error  of 

6  H 

10%  in  -r—  will  cause  an  error  of  only  1.7%  in  —  .  Measured  and 

6  'pc 

H  P 

calculated  values  of  o/pc  are  found  to  agree  within  20%  as  in  the 

P 

following  v.  le: 

Table  (4-2) 


H 

H 

fps 

— —  from  Measurements 
n  c  - 

p  of  w't'  , 

_2_  from  Eq.  (4-4), 
pcp  ft°F/sec 

ft  °F/sec 

50.0 

2.16 

2.58 

20.0 

0.95 

1.12 

10.7 

0.32 

0.312 

Equation  (4-4)  is  in  fact  based  on  the  parameters  of  the  whole  boundary 
layer,  and  may  not  take  into  account  properly  all  the  conditions  near 
the  wall . 

In  view  of  the  preceding  discussion,  it  can  be  stated  that  the 
accuracy  of  u'w'  and  w't'  measurements  is  better  than  ±  15%. 

4.3  Structure  of  the  Wall  Layer: 

In  the  previous  sections  the  results  were  presented  for  the 
boundary  layer  as  a  whole.  No  theory  treats  the  motion  in  the  whole 
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boundary  layer  at  the  same  time.  It  is  customary  to  diviae  the  flow 
field  into  regions  like  viscous  sublayer,  wall  layer,  outer  layer,  etc. 
Various  theories  reported  in  Chapter  II  have  been  developed  for  the 
notion  in  the  surface  layer  of  the  atmosphere.  The  basic  assumptions 
involved  are  those  of  plane  homogeneity  and  constant  fluxes.  It 
follows  from  the  discussion  in  section  4.1  that  plane  homogeneity  has 
been  approximately  realized  ir.  the  present  experiments,  at  least  with* 
in  6  in,  (z/S  -  0.2 i  from  the  floor.  But,  the  thickness  of  the  layer 
in  which  the  vertical  fluxes  are  approximately  constant  is  less  than 
2  in.  (z/c  '  0.06).  Experimental  data  on  velocity  profiles  in  the 
boundary  layer  and  pipe  flows  in  neutral  conditions  lave  been  shown 
previously  to  follow  the  log-law  over  a  wide  region  notwithstanding 
the  large  changes  in  the  momentum  flux.  It  has  been  pointed  out  by 
Monin  and  Yaglora  (1965)  that  the  distribution  of  average  hydrodynamic 
flow  fields  and  integral  characteristics  of  turbulent  flows  are  not  very 
sensitive  even  to  quite  significant  change  in  the  fluxes.  Some  direct 
measurements  of  the  vertical  fluxes  at  two  heights  in  the  atmospheric 
ground  layer  by  Mordukhovich  and  Tsvang  (1966)  clearly  indicate  that 
these  fluxes  change  with  height  in  the  atmosphere ,  too.  For  this  reason, 
the  data  from  a  layer  0.01  <_  -j  <  0.15,  has  been  considered  in  the 
following  discussion  of  the  structure  of  the  thermally  stratified  wall 
layer. 


4.3.1  Comparison  of  the  Mean  V e locity  and  Temp irat ure 
Profi 1 es  With  the  Similarity  Theory 

In  Figs.  31  and  32  are  presented  the  measurements  of  the 
mean  velocity  and  temperature  in  the  wind  tunnei  in  terms  of  the  coor¬ 
dinates  of  Monin  and  Obukhov's  (1954)  similarity  theory.  There  is  no 
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doubt  that  the  stability  parameter  correlates  very  well  the  uata 
pertaining  to  different  stability  conditions  in  agreement  with  F.qs. 
(2-23)  and  (2-24) .  The  choice  of  zTef  =  “^q  has  been  made  for  con¬ 
venience  so  that  zref  falls  within  the  thickness  of  the  wall  layer 
for  all  the  three  ambient-velocity  cases. 

The  temperature  data  show  more  scatter  than  the  velocity  data. 

It  can  be  seen  from  Fig.  31  that  the  log-linear  law  given  by  Eq.  (2-25) 
with  a  value  for  8  =  10  fits  the  data  best.  This  value  for  the  empir¬ 
ical  constant  8  ,  is  to  be  compared  with  the  value  of  7  reported 
by  McVehil  (1964),  8.5  by  Gurvich  (1965)  and  9.91  by  Zilitinkevich  and 
Chalikov  (1968)  from  measurements  in  the  atmospheric  surface  layer. 

From  an  earlier  set  of  measurements  in  the  same  wind  tunnel,  Plate  and 
Lin  (1966)  obtained  8=7  for  stable  case.  But,  the  values  of  u* 
used  by  them  were  determined  from  Ludwieg  and  Tillmann's  (1950)  formula 
which  gives  consistently  higher  values  than  the  measured  ones,  as  the 
flow  becomes  more  and  more  stable,  (see  Table  4-1).  A  correction  of 
u*  in  the  data  of  Plate  and  Lin  will  make  the  value  of  8  about  the 
same  as  found  in  the  present  study. 

Figure  (32)  shows  that  the  temperature  data  are  also  in  good 
agreement  with  the  log-linear  law  given  by  Eq.  (2-26).  The  curves 
representing  Eq.  (2-26)  for  8^  =  15  and  20,  are  shown  in  the  same 
figure.  It  is  seen  that  8^  =  17  will  give  the  best  fit.  That  such  a 
high  value  should  be  necessary  for  8T  is  surprising;  but  it  is  con¬ 
sistent  with  the  theoretical  conclusion  (see  Ellison  (1957)  that  the 

K 

ratio  -r^-,  which  for  the  assumptions  of  constant  fluxes  and  log- 
linearity  is  approximately  given  by 
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(4-7j 
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then  velocity  and  temperature  profiles  are  given  by 


and 


U  (z)  -  U  (  zref  ) 

A' 

Z 

P 

Zref| 

P 

u« 

px 
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- 

T  (z)  -  T  (  zre,  ) 

B' 

r 

' 

z 

T 

q 

Zref 

\q  ' 

T* 

q 

L 

• m 

(4-8) 


(4-9) 


Except  for  the  inclusion  of  an  additional  parameter  L  in  the 
power  lau  given  here,  it  is  similar  to  that  proposed  by  Deacon  (1949). 
Values  of  empirical  ">nstants  obtained  from  Figs.  35  and  36  are: 

A'  =  4.2,  p  =  0.30,  B'  =  5.7,  anc  =  0.31;  these  remain  approximately 
constant  over  a  fairly  wide  range  of  z/L  say,  0.01  z/L  ^  0.25. 

The  value  of  the  exponent  p  in  Fig.  35  is  to  be  compared  with  the 
values  obtained  in  an  earlier  study  by  Plate  and  Lin  (1966)  and  those 
observed  in  the  atmosphere  by  beacon  (1949),  who  shows  that  (1-p)  "*•  1, 
in  near  neutral  conditions,  and  is  consistently  less  than  unity  in 
stable  conditions  (e*g*>l-p  =0.73  when  IT  =0.09).  Values  of  the 
coeffecients  A'  and  B'  and  exponents  p  and  q  in  Eqs .  (4-6)  and 
(4-7),  in  fact,  vary  with  z/L  as  seen  from  Figs.  35  and  36.  The  change 
is  particularly  noticeable  for  small  z/L.  This  is  in  agreement  with 
the  results  of  Deacon. 

4.3,2  Comparison  of  the  Turbulence  Characteristics  with 
tTTe  Simi  lari  ty  Theory : 

Si 

rhe  ratio  -r. —  as  calculated  from  measured  fluxes  is 

>1 

shown  in  fig .  37a  as  a  function  of  stability  parameter 


It  is  seen 
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to  decrease  with  —  just  as  expected  from  theoretical  and  physical 

considerations  (see  e.g.,  Ellison  (1957)  and  Stewart  r1959)).-^ —  would 

S  ^ 

also  be  given  by  the  ratio  —  ,  (Fig.  37b)  if  the  momentum  and  heat 

fluxes  or  their  ratio,  remain  constant  with  height.  Since  this  is  not 

exactly  so,  the  results  of  Figs.  37a  and  37b  differ  slightly.  In  Fig. 

37b  is  also  represented  Eq.  (4-5)  which  is  based  on  the  approximation 

that  S  and  R  vary  linearly  with  -j- .  It  is  seen  that  agreement 

of  Eq.  (4-5)  with  measured  values  of  Si/  S  is  not  good.  On  the  other 

KH 

hand,  an  approximately  constant  value  of  -y—  -  0.74  corresponding  to 

Si 

constfnt  fluxes  in  Eqs.  (4-8)  and  (4-9)  is  found  within  the  range  of 

validity  of  power-law  formulae  with  constant  exponents. 

There  has  not  been  any  close  agreement  in  the  various  experimen- 
KH 

tal  determinations  of  —  ev.v.  for  near  neutral  conditions.  For 

Si 

example,  values  ranging  from  0.8  to  1.4  have  been  reported  (Monin  and 

Si 

Yaglom  (1965)).  Experimental  evidence  on  the  variation  of  -y—  with 

Si 

the  stability  is  even  more  lacking  for  stable  conditions.  Most  of  the 


estimates  have  been  made  indirectly  from  measured  wind  and  temperature 
profiles.  Direct  determination  from  the  flux  measurements  in  the 
laboratory  hav^  been  reported  by  Ellison  and  Turner  (1960)  who  conducted 
experiments  in  the  mixing  layer  of  the  salt  and  fresh  water  with  negative 

■Si 

salinity  gradients  in  a  rectangular  tube.  Values  of  - —  plotted  against 

Si 

R.  by  these  authors  show  considerable  scatter,  but  with  a  definite 

1  Si 

tendency  for  y -  to  decrease  as  stability  increases.  Fhe  present  re- 

S  Si 

suits  of  y-  aie  in  agreement  with  .hose  of  Ellison  and  Turner.  The 

M  Si  . 

initial  shaqi  change  in  tor  small  values  of  -  is  probably  due  to 

M  *■ 

the  effect  of  the  closeness  of  the  wail  as  shown  by  Fig.  30,  rather  than, 


i  < 


to  the  effect  of  stability.  Measurements  of  Mcrdukhovich  and  Tsvang 
(1966),  and  Record  and  Cramer  (1966)  in  the  near  ground  layer  of  the 

Si 

atmospher'  also  indicate  values  of  less  than  unity  for  stable 

Si 

conditions.  The  latter  are  compared  with  the  wind  tunnel  results,  in 


Fig.  37a. 


According  to  the  Monin-Obukhov  similarity  theory,  and  R^. 


must  be  universal  functions  of  —  ,  Figures  38  and  39  clear’)'  show 

2 

that  well-defined  unioue  functions  of  -r-  exist  for  R  and  Rr  in 
the  present  measurements.  Whether  these  functions  are  really  universal 
can  be  seen  only  after  comparing  data  from  different  sources.  No  such 
data  are  available  from  other  laboratory  measurements.  Some  measure¬ 
ments  in  the  atmosphere  have  been  made  by  Gurvich  as  reported  by  Monin 
and  Yaglom  (1965),  Mordukhovich  and  Tsvang  (1966),  and  Cramer  (1967). 

It  is  seen  from  Fig,  38  that  Gurvich 's  curve  agrees  with  the  present 
data  very  well.  Data  prints  due  to  Mordukhovich  and  Tsvang,  in  spite 
of  the  large  scatter,  do  not  show  any  consistent  deviations  from  what 
may  be  considered  as  a  universal  curve.  In  that  case,  it  is  immaterial 
whether  one  uses  or  R^  as  a  similarity  parameter. 

The  results  of  the  measured  turbulent  intensities  in  the  wall 


layer  are  presented  in  Fig.  40  in  terms  of  similar! ty  coordinates.  A 

'J  u  1  *”  V  V  ’ *"  ^  V 

slight  tendency  for  -  .  —  and  -  to  decrease  a. 

u,  u,  u. 

R.  increases  is  observed.  The  wind  tunnel  measurements  are  compared 
with  the  data  from  the  atmospheric  surface  layer  reported  by  Mordukhovi  h 
and  Tsvang  (1966)  in  Figs.  41,  4„’t  and  43.  Separate  symbols  have  been 
used  for  the  data  from  different  heights  to  bring  out  any  height  depen¬ 
dence  besidec  what  has  luw  n  implicit  in  K.  .  The  following  observa¬ 
tions  can  lie  made  from  this  comparison.  First,  similarity  between  the 
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wind  tunnel  data  and  the  atmospheric  data  is  noteworthy.  Second, 
atmospheric  data  show  too  much  scatter  to  give  any  precise  form  for  the 
universal  functions  f i ,  f2*  fj,  etc.,  in  Eqs .  (2-30),  (2-31),  (2-32), 
etc. _ Finally,  there  is  a  marked  height  dependence  on  the  results  of 

SP  Sp 

—  and  —  ,  both  in  the  wind  tunnel  and  the  atmosphere.  These 
u*  ** 

observations  are  generally  borne  out  also  by  the  field  measurements  at 

Massachusetts  Institute  of  Technology  reported  by  Cramer  (1967) .  The 

sf 

results  of  —  in  the  present  study  are  also  in  good  agreement  with 
an  earlier  set  of  measurements  (see  Fig.  42)  made  in  the  same  wind 
tunnel  and  at  the  same  station  by  Cennak  and  Chuang  (1965) .  They  did 
not  measure  u*  ,  however,  and  therefore  their  data  for  only  *  10  fps 
could  be  used  for  comparison  here,  using  value  of  u„  as  obtaine  1  from 
the  present  experiment.  Earlier,  Monin  (  962)  reported  a  value  of 


\yw  /u„  =■  0.7  in  the  neutral  air,  which  is  inconsistent  with  a  value 
of  about  1.3  indicated  by  the  measurements  discussed  here,  as  well  as 
vith  other  measurements  reported  by  Panofsky  and  McCormick  (1960) ,  and 
King  (1965).  Monir.’s  estimate  was  not  based  on  the  directly  rieasured 
value  of  u*  ,  which  might  be  the  reason  for  such  difference.  Some 


differences  in  the  reported  values  of  u'  /u„  are  noted  which  may 
partly  be  due  to  the  fact  that  even  under  neutrei  conditions  \J  u'^/u* 
varies  with  height  unless  the  free  stream  turbulence  level  is  very  high 
(see  Plate  and  Sandborn  (1966)). 

There  is  some  contradiction  between  the  shapes  of  the  universal 

j —— 

function  \/t'VT„  *  f.,  (  -■  ) ,  in  stable  conditions  as  reported  by 
Monin  (1962),  who  shown  that  it  decreases  with  an  increase  in 


stability, 


licated  by  *hc  wind  tunnel  and  the  atmospheric 


measurements  represented  in  fig.  43,  which  show  ihat  or  the  same 


height  \J t'2/T  *  increases  with  an  increc  ;e  in  the  stability  The 
only  other  data  due  to  Cramer  (1967)  favors  the  latter.  Clearly,  more 
dat"  :ith  accurate  heat  flux  measurements  is  needed  to  resolve  this 
serious  difference. 

On  the  whole,  it  can  be  observed  that  the  Monin-Obuknov  similar¬ 
ity  theory  is  well  supported  by  the  wind  tunnel  measurements,  and  this 
theory  should  therefore  form  the  basis  for  laboratory  modeling  of  the 
near  ground  layer  of  the  atmosphere. 

4.3.3  Comparison  of  the  Present  Measurements 
with  Ellison's  Theory: 

The  results  of  the  present  study  can  be  used  to  calculate 
the  decay  times  T^,  T2>  and  ,  using  Eqs.  (2-33)  through  (2-35), 
and  the  lengths  and  1^  introduced  by  Ellison  (1957)  as  discussed 

in  Chapter  II.  The  results  of  these  calculations  are  summarized  in 
the  following  tabl?. 
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Table  (4-3) 


z_ 

jC 

R. 

l6 

Rf 

h 

T! 

~1 

V 

Si 

Si 

0 

T2 

T3 

w2 

T  JT7 
*2 

ft 

ft 

.011 

.0033 

.34 

5.18 

5.72 

1.96 

0.55 

.048 

.030 

.025 

.0075 

.43 

5.36 

6.01 

2.57 

0.42 

.044 

.10 

.021 

.77 

5.58 

6.47 

4.97 

0.21 

.036 

.011 

.006 

'8 

5.08 

5.16 

1 .43 

1.01 

.081 

.055 

.02a 

.0’ 

.38 

5.42 

5.14 

1 .95 

0.64 

.ro5 

.10 

.033 

.67 

4.52 

5.93 

3.99 

0.30 

.051 

.011 

.0078 

.25 

4.57 

4.86 

1.23 

1  .33 

.110 

.073 

.025 

.013 

.38 

4.88 

4.88 

1.86 

0.82 

.081 

.10 

.039 

.68 

4.93 

5.77 

3  Ci 

0.38 

.062 

Ellison 

assumed  CT^/ i 2 )  * 

1,  in 

Eq.  (2-36) 

to  predict  a  value  of 

R- 

fer 

*  0.15. 

The  experimental 

results  indicate 

that  T 

j/Tj  ,  as  well 

as  q'  /w'  vary  significantly 

£ 

with  ~r  .  In  so  for  as  Rf 

6  rcr 

for  near  critical  conditions  as 

with  stability  and  also, 

is  equal  to  the  value  of 

indicated  by  Eq.  (2-36), 

to  some  extent . 

/  T,  “72\-i 

1  ♦  J- 

T,  ~1  ’ 

\  2  w' 

the  present 

data 

are  in  ag: 

reement 

with  Ellison's  assertion  that  R,. 

ter 

may  be 

even 

less 

than  0.1s 

The 

assumption 

of  the 

const  incy 

T 

of  1/T 

,  is  no 

supported  by  the  experiments,  however.  But,  —  15  seen  to  remain 

3 

fairly  constant,  ano  not  much  different  from  a  value  of  n  used  in 
Lq.  (2-37)  to  predict  a  value  of  ,  3  0.3  for  near  neutral 

conditions.  This  is  also  supported  ty  the  results  in  Fig.  44.  The 
length  scales  1^,  and  I,  promised  hy  Fllison,  decrease  with  an 
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The  noraalited  spectra  for  *  30  fps  are  shown  in  Figs.  45  arc!  46. 

The  shift  in  the  spectra  with  the  position  in  the  boundary  layer  is  quite 

significant.  As  increases,  the  proportion  of  the  energy  contained 

in  high  wave  numbers  becomes  less  and  less .  The  high  wave  number  end 

_  *? 

of  all  the  spectra  shows  „  k  behavior  in  agreement  with  Heisenberg’s 

/x 

( 1948}  theory.  An  extensive  subrange  showing  k  '  behavior  is 
observed  for  U  *  30  fps,  when  s/5  is  large.  The  exten*-  of  this 
subrange  becomes  less  and  less  as  U  or  -/6  is  decreased.  As  has 
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been  discussed  by  Alkseev  and  Yaglos  (1967),  determining  the  extent  of 
inertial  subrange  froa  the  measured  one -dimensional  spectrum  is  mis¬ 
leading  in  the  sense  that  since  a  one -dimensional  spectrun  varies  more 
smoothly  than  the  corresponding  three-dimensional  spectrum,  the  former 
may  be  approximated  by  the  formula  F  (kj)  k^”^^  ,  over  a  consider¬ 
ably  greater  range  of  wave  numbers  than  the  three-dimensional  spectrua 
E  (  k  ). 


4.4.2  Comparison  Kith  Kolmogorov's  Similarity  Theoi 


In  order  to  co:cpare  the  observed  spectra  with  Kolmogorov's 


(1941)  similarity  theory,  it  is  necessary  to  express  them  in  the 

S  1/4 

reference  frame  of  similarity  coordinates,  i.e.,  F  (k.)/(ev  )  ' 
k  1/4  3/4 

versus  l/kg  where,  k  =  e  /v  ,  is  the  Kolmogorov  wave  number. 

Commonly,  e  is  determined  from  the  u’ -  spectrum  by  using  che  isotropic 


relation 


e 


(kj)  dkx 


(4-12) 


This  method  is  not  available  in  the  present  study  since  u'-  spectrum 
could  not  be  measured  because  of  difficulties  in  obtaining  a  signal 
proportional  to  u'  alone.  Even  otherwise,  it  is  doubtful  whether 
Eq.  (4-12)  will  be  valid  for  sufficiently  stable  conditions,  where 
the  effect  of  low  velocities  as  well  as  buoyancy,  is  to  make  a  larger 
part  of  the  spectrum  deviate  from  the  assumptions  of  local  isotropy. 

As  discussed  in  Chapter  II,  for  the  region  in  which  diffusion 
terms  in  the  energy  equation  will  be  negligible,  Eq.  (2-18)  can  be  used 
for  calculating  c  .  Accordingly,  e  has  been  determined  (see 


Table  VI)  for  the  layer  .009  <_ —  <_  .075  and  the  nondisens ionali zed 

spectra  Fyt  {  —  )  *  Fyt  (kj)  /  (evs)1/f4  .and  Fw,  (  )  = 

s  s 

S  1/4 

F  ,  (k  )  /  (  ev  )  plotted  in  Figs.  <7  and  48.  It  is  seen  that  the 
spectra  pertaining  to  different  velocities  (stabilities)  and  to  differ¬ 
ent  positions  in  the  boundary  layer  fall  close  together  in  the  high 
wave  number  region,  and  branch  off  different  points,  depending  upon 
the  size  of  energy  containing  eddies,  from  what  may  be  considered  a 
universal  curve.  There  is  also  ooserved  a  small  but  consistent  shift 
in  the  spectra  with  z/6  ,  such  that  the  constant  o  in  the  spectral 
form  for  the  inertial  subrange  given  by 


(4-13) 


increases  with  an  increase  in  z/6.  This  shift  cannot  be  explained 
with  certainty.  If  the  spectra  for  30  and  20  fps  are  compared  at  the 
same  z/6  ,  no  significant  change  in  b  is  noticeable.  Therefore,  the 
observed  shift  is  not  due  to  the  effect  of  stratification,  but  perhaps 
due  to  the  error  in  the  estimation  of  e  .  A  more  accurate  determina¬ 
tion  of  e  will  be  needed  to  check  whether  there  is  really  some  effect 
of  stability  on  the  Kolmogorov  constant. 

The  present  experimental  results  can  be  compared  with  the  spec¬ 
tral  form  obtained  from  Heisenberg's  (1948)  theory  which  for  three- 
dimensional  energy  spectrum  in  the  equilibria  range  can  be  written  as 
(see  Rotta  (1962)) 


E  (lr)  = 


-4/3 


(4-14) 


t>4 


2/3  k-5/3 

f  c  3  .4  ' 

1  ♦  8v,  k 

l9 

T  2 

3o^  e 

in  which  is  a  constant.  After  substituting  in  Eq.  (4-14)  the 
numerical  value  of  corresponding  to  the  value  of  0.48  fcr  the 
Kolmogorov  constant  a^  for  u'  spectrum,  and  transforming  it  in 
terms  of  universal  parameters,  one  obtains  from  Eq.  (4-14) 


E  (n)  =  1.466 


-5/3 


1  +  10.67  n 


-4/3 


(4-15) 


in  which,  n  =  y  .  The  corresponding  one-dimensional  spectrum  for 
s 

lateral  or  vertical  velocity  fluctuations  can  be  obtained  by  using  the 
isotropic  relation  (see  Batchelor  (1953a)) 
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which  in  terms  of  n  and  =  y-  ,  can  be  written  as 
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k  k 

where  F,  (n, )  =  F  ,  ( -r^-  )  =  F  ,  ( ) 
21  v'  k  w'  v  k  ’ 
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After  substituting  from  Eq.  *' 4 - 1 5 )  into  Eq.  (4-17),  the  following 
equation  is  obtained. 
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Eq.  (4-18)  represents  Heisenberg's  one -dimensional  energy  spectrum  of 
lateral  or  vertical  fluctuations.  A  numerical  integration  of  Eq.  (4-18) 
was  performed  which  has  been  represented  in  Figs.  47  and  48  for  com¬ 
parison  with  the  experimental  data.  It  can  be  seen  that  Eq.  (4-18) 
describes  the  measured  spectra  in  the  equilibrium  range  quite  well. 


4.4.3  Spectra  of  Temperature  Fluctuations: 

The  measured  one -dimensional  spectra  of  temperature 
fluctuations  are  shown  in  Figs.  49  and  50.  These  have  been  normalized 
such  that 


J  V  (ki>  dki  ■ 1 


(4-19) 


The  same  trend  with  z/6  is  observed  for  t’-  spectra  as  for  v'-  and 
w'-  spectra.  The  high  wave  number  end  shows  k  ^  behavior  in  agree¬ 
ment  with  Corrsin's  (1951)  prediction.  The  inertial  subrange  region  is 
almost  absent  for  the  case  of  ■  10  fps;  it  is  significant,  however, 

for  U  =>  20  fps  and  z/f  -  0.37. 

00  * 
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In  none  of  the  v'  ,  w'  ,  and  t'  spectra,  could  the  so-called 
buoyancy  subrange  be  identified.  This  was  expected,  of  course,  from 
I.umley's  (1965a)  discussion  according  to  which  the  byoyancy  subrange 
would  not  occur  except  in  very  strong  inversions  where  would  be 
close  to  unity. 

The  spectrum  of  temperature  fluctuations  can  be  used  to  calculate 
e^.,  by  using  the  isotropic  relation  (see  Hinze  (1959)). 


et,  =  6  a 


♦t,  (>9 


(4-20) 


Values  of  calculated  from  Eq.  (4-20)  are  compared  in  the  follow¬ 

ing  Table  with  those  given  by  Eq.  (2-19)  using  heat  flux  measurements. 

Table  (4-4) 


z 

A 

Et,  from 

et,  from 

O 
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65 

28.2 
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The  agreement  is  not  very  good.  The  possible  reasons  for  this  are: 
First,  an  error  in  calculating  ,  ,  from  spectral  measurements  using 
isotropic  relation  (4-20),  while  the  low  wave  number  end  of  these 
spectra  up  to  a  wave  number  as  high  as  70  is  not  even  locally  isotropic; 
second,  an  error  in  the  value  of  e^,  obtained  from  Eq.  (2-19)  due  to 
neglected  diffusion  terms.  Over  what  region  of  the  wall  layer  Eqs . 
(2-18)  and  (2-19)  can  be  used  co  determine  t  and  e^,  cannot  be 
ascertained  from  the  present  measurements. 
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Chapter  V 
CONCLUSIONS 

The  results  of  measurements  in  a  stably  stratified  turbulent 

boundary  layer  in  the  wind  tunnel  have  been  discussed  in  Chapter  IV. 

A  range  of  gradient  Richardson  numbers  from  0  to  0.25  in  th<  boundary 

layer  has  been  covered.  With  the  difference  in  the  wall  temperature 

and  ambient  air  temperature  almost  fixed  in  the  experiment,  has 

become  necessarily  correlated  with  ambient  velocity  and,  therefore ,  with 

R  .  For  this  reason,  the  Reynolds  number  effect,  though  relatively 
*6 

small,  also  has  shown  up  when  the  effect  of  change  in  on  turbulent 
quantities  is  studied.  The  rms  fluctuations  and  turbulent  fluxes  as 
normalized  by  and  (T^-T  )  are  shown  to  decrease  considerably  with 
increase  in  stability.  The  vertical  component  is  relatively  more 
affected  than  the  longitudinal  component. 

Special  attention  has  been  given  to  the  lowest  (excluding  the 
viscous  sublayer)  15%  of  the  boundary  layer,  i.c.  the  wall  layer.  The 
results  obtained  for  this  layer  have  been  compared  with  Monin  and 
Obukhov's  similarity  theory  and  other  theoretical  models.  The  theoret¬ 
ical  requirement  of  constant  fluxes  is  met  at  roost  by  only  a  third  of 
the  wall  layer,  but  plane-homogeneity  in  the  sense  that  the  flow  char¬ 
acteristics  do  not  change  noticeably  in  the  direction  of  flow  has  been 
realized  at  the  measuring  station. 

Mean  velocity  and  temperature  profiles  in  the  wall  layer  are  in 
fair  agreement  with  Monin  and  Obukhov's  (1954)  similarity  theory. 
Logarithmic- linear  laws  with  K  *  10  and  fci.j.  *  17  in  Fqs .  (2-25)  and 
(2-2 6)  represent  the  mean  velocity  and  temperature  data  quite  well. 
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Significant  departures  from  the  corresponding  linear  relations  for 
kz  3  U  z  3T 

S  = —  t—  ,  and  R  ,  however,  are  noticed.  These  are  shown 

U*  dZ  T*  dz  * 

to  be  better  represented  by  power-law  relations  which  also  suggest  the 
power-law  form  of  the  velocity  and  temperature  profiles  as  in  Eqs.  (4-8) 
and  (4-9) . 

KH  2 

The  parameters  ,  R.  and  R.  are  unique  functions  of  -r- 

*M  1  *  1 

in  agreement  with  the  similarity  theory.  Close  agreement  with  the 

limited  measurements  that  are  available  from  the  atmospheric  surface 

layer  further  shows  that  these  functions  may  well  be  universal.  Since 

R.  is  a  universal  function  of  -r-  ,  it  is  equally  valid  to  use  the 

1  L 

former,  which  is  more  easily  determined,  in  the  similarity  representation 

z  Si 

in  place  of  — .  The  ratio  is  consistently  less  than  unity  (up  to 

0.6  in  the  experiment)  in  stable  conditions  so  that  the  often  made 

assumption  of  -rr-  *  1  may  not  be  justified  in  even  moderately  stable 
*M 

stratification.  _  _  _ 

Turbulent  quantities  - , - and  -  are  slowly 


--  and 


varying  functions  of  R^  .their  magnitude  decreasing  with  increasing 

stability.  In  magnitude  these  compare  well  with  atmospheric  data  which 

show  the  trend  with _ R.  less  clearly  because  of  large  data  scatter. 

The  quantity  —  ,  although  somewhat  dependent  on  height,  is  seen 

*  * 

to  increase  with  increasing  R^  .  This  agrees  with  atmospheric  measure¬ 
ments  reported  by  Cramer  (1967),  but  disagrees  with  those  of  Monin 

(1962)  who  shows  -= —  as  a  decreasing  function  of  ~  .  In  the  latter, 

*  L 

however,  T,  was  not  determined  from  direct  heat  flux  measurements,  and 
it  is  possible  that  this  resulted  in  such  p  difference.  More  data  is 
needed  to  resolve  th;s  serious  difference. 
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In  general,  it  can  be  concluded  that  mean  flow  and  turbulent 
characteristics  of  the  wall  layer  in  the  wind  tunnel  are  well  represent¬ 
ed  by  the  Monin-Obukhov  similarity  theory  and,  therefore,  this  theory 
should  form  a  basis  for  modeling  of  similar  characteristics  of  the 
atmospheric  surface  layer.  Moreover,  wind  tunnel  data  can  give  the 
form  of  various  universal  functions  in  the  similarity  equations  more 
precisely  than  the  atmospheric  data  can  do  because  of  large  scatter 
and  variable  conditions  in  the  latter. 

The  results  of  Ellison's  (1957)  theory  as  to  the  general  trend 


_KM  ’ 


of  variation  of  —  *  Yw't'  *  Hi  ant*  St  stability  aii  verified 


by  the  present  experiment.  The  assumption  of  the  constancy  of  the 

T  "72 

ratios  and  3-y  in  his  theory,  is  not  borne  out  by  the  measure- 

9  u  *  t*  — — rx 


4 

ments  which  show  both  =—  and 


S2. 

7-2 


increasing  with  increase  in 


stability.  The  assumed  numerical  value  of 


1  Q 


,2 


T_  tT*2 


=  5.5,  however,  is 


in  agreement  with  the  experimental  trend  for  strong  stability  conditions, 

and  consequently,  Ellison's  prediction  of  R^.  ^0.15  appears  to  be 

cr 


roughly  correct. 

The  normalized  roectra  of  velocity  and  temperature  fluctuations 
show  strong  dependence  on  —  and  .  From  the  non-dimensionalized 
representation  of  the  spectra  in  terms  of  Kolmogorov's  similarity 
theory  coordinates,  it  can  be  concluded  that  all  spectra  have  a  shape 
which,  irrespective  of  thermal  stratification,  follows  a  universal  curve 
for  large  wave  numbers  that  correspond  to  the  equilibrium  range.  The 
magnitude  of  the  energy  in  this  range  itself  gets  reduced,  however,  due 
to  the  combined  effect  of  decrease  in  Reynold  number  and  increase  in 
stability.  Stratification,  together  with  shear,  affects  the  spectra 
also  m  the  range  of  energy  containing  eddies.  Hut  no  particular 
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subrange  is  indicated  in  which  buoyancy  will  act  to  the  exclusion  of 
shear  to  give  the  so-called  buo>ancy  subrange.  This  is  in  agreement 
with  Lumley's  (1965a)  discussion  of  the  conditions  for  the  existence 
of  a  buoyancy  subrange.  Although  the  small  scale  structure  of  tur¬ 
bulent  motion  remains  unaffected  by  stability,  the  overall  energy  level, 
and  magnitude  of  turbulent  production  and  dissipation  get  considerably 
reduced  as  stabiltiy  increases. 

In  summary,  it  can  be  concluded  that 

1.  The  effect  of  increasing  stability  is,  in  general,  to  sup¬ 
press  turbulence  in  the  boundary  layer  and  to  make  the  flow 
more  anisotropic. 

2.  The  structure  of  the  wall  layer  is  well  represented  by 
the  Monin-Obukhov  similarity  theory.  In  particular,  the 
logarithmic- linear  law  for  mean  velocity  and  temperature 
profiles  with  the  value  of  empirical  constants  about  the 
same  a?  has  been  found  from  atmospheric  measurements,  is 
valid.  It  is  also  shown  that  a  power-law  similarity  profile 
can  give  even  a  better  fit. 


H 

3.  In  the  wall  layer,  parameters  ,  fT  and  R^.  are 

universal  functions  of  the  stabUity  ratio  ~  ,  and  so  are _ 

JP"  -T^  A"?' 

the  turbulent  quantities  ’ - ,  - ,  - ,  -= — , 

M  u*  ’>*  u*  T* 

etc.,  in  agreement  with  the  results  of  the  similarity  theory. 
Kil 

—  drops  sharply  away  from  the  wall  and  is  consistently  less 

Si  ,r— T 

than  unity  in  stable  stratification.  V  u'"/u,  and 


V  t’VT*  also  show  some  height  dependence  over  and  above 


that  considered  in  the  stability  parameter. 
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4.  Good  agreement  of  the  present  laboratory  data  for  stable 
conditions  with  the  Monin-Obukhov  similarity  theory,  which 
has  been  shown  previously  by  others  to  represent  the  struc¬ 
ture  of  the  atmospheric  surface  layer  shows  that  this  theory 
should  form  the  basis  for  simulation  of  near  ground  atmos¬ 
pheric  motion  by  wind  tunnel  flows. 

5.  Kolmogorov's  local  isotropy  theory  is  valid  for  turbulent 
spectra  in  the  equilibrium  wave  number  range  irrespective 
of  stratification.  The  extent  of  this  range,  however,  is 
greatly  reduced  as  stability  increases. 

6.  No  byoyancy  subrange  is  indicated  by  the  measured  spectra 
of  V '  ,  w 1  and  t ' . 

The  following  recommendations  can  be  made  for  any  future  investi¬ 
gation  of  the  thermally  stratified  boundary  layer  from  the  point  of  view 
of  atmospheric  simulation. 

1.  Measurements  of  the  mean  flow  and  turbulent  quantities  in 
unstable  stratification,  as  well  as  more  strongly  stable 
conditions  than  covered  by  the  present  experiment,  should  be 
made  to  study  whether  conclusions  of  the  present  study  can 
be  extended  to  include  all  the  stability  range  of  interest. 

2.  Efforts  should  be  made  to  thicken  the  constant  flux  layer 
in  the  wind  tunnel  by  the  use  of  surface  roughness  or  shear 
grid,  to  make  the  flow  conditions  more  near  to  those  of  the 
theoretical  model. 

Attempts  should  be  made  to  make  R.  as  independent  of 
R  as  practicable  by  varying  the  temperature  or  the  wall 
and  that  of  ambient  air. 


3. 
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4.  The  dissipation  e  ,  should  be  determined  by  using  a  more 
direct  method  e.g.,  from  the  spectrum  of  longitudinal 
fluctuations.  The  fact  that  the  normal  wire  is  sensitive  to 
to  both  velocity  and  temperature  fluctuations  make  this 
somewhat  difficult  and  new  experimental  techniques  need  be 
developed  for  this  purpose. 
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APPENDIX  A 
YAWED  WIRE  RESPONSE 

For  a  given  hot-wire  anemometer  set  up,  the  voltage  output  across 
the  wire  must  depend  on  the  total  velocity  Up  ,  angle  of  yaw  6  , 
and  the  difference  in  the  temperature  of  the  wire  and  that  of  local 
fluid.  So  that  in  functional  form  one  can  write 

E  -  f  (UTot,  6  ,  T)  .  (A-l) 

For  convenience,  two  variables  U_  .  and  6  can  be  combined  to 

Tot 

constitute  what  may  be  called  the  'effective  velocity,'  »  for 

heat  transfer.  So  that 


E  -  f  (Uef£  .  T) 

and  in  differential  form 


CA-2) 


Recognizing  the  fact  that  in  addition  to  the  component  of  the 
total  velocity  which  is  normal  to  the  wire,  the  one  parallel  to  it  also, 
affects  the  heat  transfer  from  the  finite  wire,  Dime  (1959)  and  later 
Webster  (1962)  suggested  the  following  expression  for  in  the 

absence  of  turbulence. 


a 


C05  ‘  ') 


( A  -  4 ) 
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In  Eq.  (A-4) ,  a  is  an  empirical  constant  having  a  value  between 
0.1  and  0.3.  Webster  determined  an  average  value  of  a  =  0.2.  His 
measurements  show  considerable  scatter  with  values  of  a  ranging  be¬ 
tween  0.1  and  0.3,  but  no  systematic  variation  with  length-to- 
diameter  ratio  of  the  wire.  Precise  measurements  of  heat  transfer  from 
hot-wires  were  also  made  by  Champagne,  Sleicher  and  Wehrn  in  (1967) 
who  showed  that  Eq.  (A-4)  correlates  the  data  very  well.  From  this 

study  value  of  a  is  found  to  depend  primarily  on  l'/d  ratio.  For 

1 1 

platimum  wires,  a  is  approximately  0.20  for  -g-  =  200,  decreases 

1 1  i  » 

with  increasing  -j-  and  becomes  effectively  zero  at  -j-  =  600.  No 
significant  differences  in  the  values  of  a  are  found  for  wires  of 
different  materials. 

In  the  turbulent  field,  referring  to  Fig.  5,  one  has  in  place  of 
Eq.  (A-4) 

U4ff  =  |(U  ♦  u')4  *  w’2}  sin2  (6  ♦  dfi)  ♦ 

a4  | (U  ♦  u')2  +  w'2 

Using  the  trignometrical  identities  for  sin  (0  ♦  dO)  and 

cos  (0  ♦  dO)  ,  and  further  recognizing  that 


cos 


(6 


de) 


„2 


(A-5) 


sir.  do 


COS  d1' 


w ' _ 

((U  ♦  u’)2  ♦  w'2}  1/2 

i  • 

U  ♦  u ' 

I  ■* 

s ,  (  U  ♦  t!  ‘  )  *  » 

i 
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the  following  equation  can  he  obtained  after  neglecting  second  order 

„  u'  v'  ,  w' 

terms  m  -jj- ,  -jj— ,  and  —  : 

(1  ♦  a2  cot2  6)  (1  ♦  2  y-)  ♦  2  (1  “  a2) 

(A-6) 

Equation  (A-6)  after  differentiation  gives 


Ueff  *  11  s*n  0 


IT  COt  e 


1/2 


dU 


eff 


=  (1  ♦  a  cot  6) 


-1/2 


,  +  2  —  >  2  f  — 1  Cl  -  »  3  ..cot  e. 

U  (  «  3  (1  ♦  a2cot2e) 


-,-1/2 


(1  ♦  a2  cot2  6)  du’  sin  6  (1  -  a2)  dw 


'  cos  6 


(A-7) 


After  Bemaulli  expansion  of  the  second  parenthesis  terms  and 

u*  w 1 

again  neglecting  second  and  higher  order  terms  in  —  and  —  and  also 
terms  like  u'du*  and  u'dw'  ,  one  obtains 


dU^^  *  ( i  ♦  a2  cot2  0)  ^2 


2  >  2 
(1  ♦  a  cot“  0)  du'  sin  <:♦(!-*)  dw'  cos  0 


(A-8) 


After  substituting  fro*  tq.  (A-8)  into  F.q.  (A-3)  and  making  use 
of  the  common  assumption  that  for  saaJl  fluctuations  differentials  can 
be  replaced  by  fluctuations  themselves,  the  following  equation  is 


obtained 
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,  3E  ,  .  2  .  2  2  .,1/2 

e*  =  -rrr  (sin  0  +  a  COS  0) 

eff 


...  .  (1  -  a  )  cot  0 

u  +  - 7 - 2 — 

(1  +  a  cot  0) 


♦  f 

3T 


(A-9) 


Since  the  wire  sensitivity  is  determined  from  calibration  in  a 
turbulent  free  stream,  it  follows  from  Eq.  (A-4)  that 


3E  r  ■  2  n  2  2  ..  3E 

3U  „  (S111  9  +  3  COS  9)  =  gjj 

eff 


(A-10) 


so  that 


e '  * 


3E 

3U 


u'  ♦ 


(1  +  a2  cot2  0) 


w'  cot  9 


♦  t- 
3T 


(A-ll) 


Equation  (A-ll)  is  the  basic  response  equation  of  the  yawed  wire, 

Comparing  it  with  Eq.  (3-1)  yields  the  required  relationship  between 

_  ....  3E  ,  3E 

the  two  sensitivities  and  ,  viz., 


L  iL  ,  (1  -  a*)  cot  9 


U  36 


(1  ♦  a‘  cot2  9) 


3E  3F. 

W  *  c  cot  9  aU 


(A- 12) 


where , 


(l  -  a") 

- , - *  - 


(1  *  a*  cot  •>] 


c  * 


(A- 13) 
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It  can  be  seen  that  c  decreases  with  decreasing  yaw  angle  which  is 
what  may  be  expected. 

Thus,  for  a  known  a  and  6  ,  c  can  be  calculated  from  Eq. 

(A- 13) .  The  following  table  gives  values  of  c  for  different  values 
of  9  ,  and  for  a  =  0.20 

6  60°  50°  45°  40°  30° 

c  0.948  0.934  0.923  0.908  0.857 

A  more  rigorous  derivation  of  the  hot-wire  response  equations  by 
Champagne  et.  al . ,  (i 967T  leads  to  a  value  of  c  =  0.925  for  0  =  45°. 
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APPENDIX  B 

COMPUTER  PROGRAM  FOR  LEAST  SQUARE  FIT 
OF  EQUATION  Y  =  A  X2  +  B  X 


•FORTRAN 


program  lstso 


ooono? 

DIMENSION  SUm (5) t  HE 0(8) » X { 12) » Y ( 12) 

000002 

FlNIS**6HFIMtSH 

000004 

1 

RE AO  102.HED 

000011 

!F(HE0(1)-FIMSK)  2*3q*2 

000013 

2 

PRINT  104*  H»n 

000020 

Pr  AD  100  »NP 

000025 

DC  4  1*1.5 

000027 

4 

SUM ( 1 ) >0  * 

000033 

DO  10  !«1«NP 

000035 

READ101 i X ( 1 ) » Y  ( I) 

0C0044 

SUM (4)  •  SUM (4)  ♦  Y ( I )  *X(I) 

000050 

SUM  (5)  ■  SUm(5)4Y(!)*X(I)*X(I) 

000053 

00  10  J*1 *3 

000034 

10 

SUM(J)«SUM(J)  ♦  X(I)4«(J*1) 

000070 

0«SUM (3) *SUM ( 1 )  -  SUM(2)*SUM(2) 

000073 

PRINT109, (SUM(J) • J" 1 • 5 ) *D 

000106 

A« (SUM (1 ) *SUM (5) -SUM (2) *SUM (4) ) /D 

000U2 

B* (SUM (3)  *SUM<4)  •  SUM (2)  *SUM(5))/D 

OQ0116 

WRITE  (6.105)  A.B 

000125 

WRITE  (6.106) 

000130 

DIFF2»0. 

000131 

DO  20  Irl'.NP 

0Q0133 

YC"A*X(I)*X(n  ♦  3*X  ( I ) 

000137 

OIFF  ■  Y(I)-VC 

000141 

nIFF2«DlFF2  ♦  OIFF*OIFF 

000143 

20 

WRITE  (6.107)  X(I) .Y(I) .YCtDIFF 

000161 

WRITE (6* 108)  DIFF2 

000166 

GO  TO  1 

000167 

30 

call  exit 

000170 

100 

FORMAT (14) 

000170 

101 

FORMAT(2F10.O) 

000170 

102 

FORMAT ( 8 A 1 0 ) 

000170 

104 

FORMAT(lHl»8AlO) 

000170 

105 

FORMAT (3H0A«,E16.7»6X.2HB«»E16, 7) 

000170 

106 

FORMAT ( 1H0 »7x » 1HX* 1 IX . THY  GlVEN*6X . 1 OHY  COMPUTED* 
15X.10HOIFFERENCE) 

000170 

107 

F0RM*T(1h  4  E15.4) 

000170 

108 

FORMAT (30HOSUM  OF  SQUARED  DIFFERENCES  ■» 

1F16.7) 

000170 

109 

FORMAT ( lK  .6F15.4) 

000170 

END 
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TABI.E  1 


SUMMARY  OE  THE  BOUNDARY- LAYER  PARAMETERS 


m 

u„ 

fps 

B 

r~ 

°F 

6 

ft 

R  K 

EM 

u, 

fps 

T 

°F 

L 

ft 

80.  05 

n 

39.  0 

2.  23 

3.  u2xi  o;? 

.  0113 

M9H 

5,  38 

17  6 

B 1 

30.  00 

HU 

40.  0 

2.  23 

3,  02x10" 

.  0112 

IH>9 

5  27 

17  5 

78 

30.  05 

BEBrH 

2.  23 

3.  02x10* 

010 

H^ 

5  00 

210 

78 

30  40 

Wmm 

42.  0 

2.  23 

3.  05x10* 

0116 

5  28 

19  3 

78 

30  10 

11119 

40  0 

2.  23 

3.  02x10? 

0111 

W  r?M 

5  37 

18  5 

70 

30  10 

119.  2 

40.  ! 

2  23 

3.  03x10* 

0110 

5  45 

18  6 

78 

20.  25 

118.  0 

40.  0 

2.  29 

2.  09x10? 

0234 

4  19 

8  20 

78 

19  95 

116  2 

40.  0 

2  29 

2  06x10? 

0247 

0  584 

3  97 

7  95 

78 

20.  15 

111.0 

37  0 

2  29 

2  09x10" 

0236 

0  581 

3  82 

8  11 

78 

20  20 

117  5 

40  0 

2.  29 

2  09x10? 

0234 

0  568 

3  89 

7  67 

78 

19  80 

118  4 

39  0 

2.29 

2  04x10? 

0262 

0  604 

4  32 

7  80 

70 

19  80 

118  4 

39  0 

2  29 

2.  34x10* 

0262 

0  604 

4  32 

7  30 

78 

9  70 

117  0 

38  0 

2  37 

'  05x10? 

113 

0  221 

3  28 

1  40 

78 

9  35 

117  0 

38  0 

2  37 

1  01x10 

121 

0  228 

3  51 

1  40 

78 

10  10 

117  0 

38  0 

2  37 

1.  09x10? 

104 

0  241 

3  43 

1  59 

78 

10  05 

114  0 

40.  3 

?.  37 

1  08x10? 

.  0981 

0.  243 

3  26 

1  77 

76 

10  70 

_ 

115.  5 

40.  0 

2.  37 

1.  1 5x10^ 

0883 

. .  ^  .  —  i 

0.  246 

3  17 

1  81 

90 


I  AHI.I-:  I! 


Ml-: AN  VKI-OOTY  AND  TKMI'KHATI'Hi:  DATA 


Kun  No. 

*  -  78  ft 

Ut  ’  30.05 

To  »  U‘).9 
T0  -  39.0 

3  U 

in.  fps 

fps 

r 

°¥ 

Itmi  No  2 

x  -  78  ft 

II  »  30.00  fps 

T  =  119.8  °F 

T0  -  40.0  °F 

z  U  T 

in.  fps  °F 

Klin  No  3 

x  ■  78  ft 

II  '■*  30. OS  fps 

T  »  113.3  °F 

T0  ■  42.0  «>F 

i  U  T 

in.  fps  °F 

0  19 

IS  14 

91.  83 

0  19 

14.  75 

91  70 

0  125 

15  50 

84.  30 

0  34 

16  72 

94  50 

0.  25 

15  23 

92  34 

0  200 

17  36 

86.  7 

0  54 

18  14 

97.  74 

0  41 

16  34 

96  04 

0  355 

18  25 

90  0 

0.  96 

19.  35 

101.  78 

0.  66 

18.  03 

97.  65 

0  66 

19.  57 

94  8 

1  41 

20.  50 

102  70 

1.  05 

19.  33 

100.  95 

0.  94 

20.  35 

97  3 

1.  95 

21  40 

104  87 

l  50 

20.  37 

103  30 

1.  45 

21  33 

98.  5 

2  50 

21. 92 

106  43 

2.  00 

20.  S3 

104.  42 

1. 98 

22.  00 

100  5 

3.  01 

22.  43 

108.  60 

2.  60 

21.  50 

105.  70 

2.  47 

22.  58 

101.  5 

3.  47 

22.  90 

109.  13 

3.  16 

22.  30 

108.  25 

2.  98 

23.  10 

103.  0 

4.  01 

23.  23 

110.  60 

3.  88 

23.25 

109.  60 

3.  77 

23.  80 

104.  5 

4.  49 

23.  65 

111  60 

4.  70 

23.  95 

111. 25 

5.  07 

24.  67 

108  3 

5.  03 

24.  06 

111.  74 

5.  70 

24.  55 

111.  88 

6.  39 

25.  40 

5.  16 

112.  91 

6.  62 

24.  76 

2.  66 

7.  75 

26.  10 

108.  5 

6.  00 

24.  75 

113. 79 

7.  63 

25.  40 

13.  79 

9.  40 

26.  63 

109.  3 

6.95 

25.  30 

114.  61 

8.  75 

25.  85 

114. 70 

10.  63 

27.  08 

110.0 

7.  98 

25.  85 

115.48 

9.  80 

26.  22 

115.  47 

12.  00 

27.  32 

110.6 

8.  90 

26.  20 

116.  08 

10.  85 

26.  63 

116.  30 

13.  37 

27.  62 

110.  8 

9.  90 

26.  55 

116.  52 

11.93 

27.  20 

116.  87 

14.  73 

27.  93 

1112 

10.  80 

26.  83 

117.  00 

12.  95 

27.  38 

117.  38 

16.  14 

28.  37 

111.  5 

11.90 

27.  22 

117.  57 

14.02 

27.  72 

117.70 

17.  57 

28.  65 

111.  7 

12.  77 

27.45 

118.  00 

15.  05 

28.  05 

118.  00 

•.8.  82 

29.  07 

111  9 

13.  80 

27.  80 

118.43 

16.  10 

28.  36 

118.  21 

20.  20 

29.  27 

112  2 

14.  75 

28.  10 

118.  70 

17.  23 

28.  64 

118.  42 

21.  30 

29.  57 

112.  6 

15.  76 

28.40 

119.  00 

18.  30 

28.  95 

118.  66 

22.  77 

29.  83 

113  08 

16.  85 

28.  65 

119.  26 

19.  30 

29.  18 

118.  83 

24  00 

29  98 

113  17 

17.  87 

28  92 

119.  30 

20.  35 

29.  36 

119  00 

18.  90 

29.  10 

119.  30 

21  45 

29.  60 

119. 30 

19.  87 

29  23 

119.  39 

22.55 

29.  81 

119.  55 

20  55 

29  47 

119  52 

23.50 

29.  93 

119.  75 

21  60 

29.  70 

119  65 

22  65 

29.  84 

119.  74 

23  50 

29  95 

119  87 

I  \  I : I  I  II 


Ml  \\MI<m-|I\  \  I  *  i  I :  IM  I'l .  I J  A  I  i  It  I  .  DATA  Continued 


1 \. 

»  1 

j  1 t i  iii  N i ♦  ‘  i 

Run  No  6 

\  -  'S 

ft 

X 

/s  II 

X 

-  70  ft 

II  ill 

.  Ill  l|*N 

|i 

ill  III  I'ps 

>  .30.  10 

fps 

1.  -  1  I9.i.  "1 

1 

II '.Ml  °l 

T. 

*  113.2 

°F 

To  ■  «■ 

,0  1>I 

'o  -  ■ 

1(1.0  °l- 

T 

0 

*  40.0  °F 

z 

U 

2 

1 

°f 

z 

11 

T 

i 

IJ 

T 

in. 

fps 

in. 

in. 

fps 

°F  ! 

1  in- 

fps 

°F 

0.  13 

15.  08 

0.  13 

91.  50 

0  082 

14.  22 

85.  17 

0. 

144 

16.  02 

88 

42 

0.  21 

16.  19 

0.  17 

93.  00 

0.  094 

15.  23 

86.  38 

0. 

156 

16.  45 

n  i 

25 

0.  29 

17.  04 

0.  28 

94.  50 

0.  199 

16.  63 

90.  17 

0. 

340 

18.  00 

94. 

30 

0.  47 

18.  03 

0.  45 

96.  87 

0.  484 

18.  37 

93.  75 

0. 

74 

19.  55 

90. 

00 

0.  67 

19.  05 

0.  68 

98.  70 

0.  89 

19.  72 

96.  66 

1. 

13 

20.  64 

101. 

00 

1.  03 

20.  32 

1.  01 

101.  00 

1.  19 

20.  48 

98.17 

1. 

56 

21.  37 

103. 

13 

1. 43 

21.  26 

1.  39 

102.  83 

1.  53 

21.  23 

100.  25 

2. 

04 

22.  11 

104 

79 

2.  07 

22.  32 

!.  87 

104.  75 

1.  95 

21.  63 

101.  68 

2. 

54 

22.  44 

106 

34 

2.  81 

22.  98 

2.  55 

107.  42 

2.  35 

22.  27 

103.79 

2. 

96 

23.  08 

107. 

50 

3.  56 

23.  60 

3.  45 

109. 25 

2.  76 

22.  56 

104.00 

3. 

44 

23,  36 

108. 

48 

4.  35 

24.  08 

4.  42 

110.  38 

3.  13 

22,  93 

105.08 

4. 

02 

23.  88 

109. 

54 

5.  04 

24.  51 

5.  00 

112. 08 

3.  53 

23.  27 

105.  70 

4, 

59 

24.  24 

110. 

38 

5.  90 

5.  80 

113.  00 

4.04 

23.  66 

106.  62 

5. 

31 

24.  87 

111. 

48 

6.  60 

25.  42 

6.  60 

111.  08 

4.  78 

23.  97 

107.  70 

6, 

30 

25  38 

112. 

17 

7.  60 

25.  87 

7.  70 

114.  87 

5.  40 

24.  45 

109.21 

7. 

25 

25.  90 

113 

30 

8.  65 

26.  31 

8.  70 

115.  87 

5.  95 

24.  95 

109.  62 

8 

40 

26  52 

114 

00 

9.  70 

26  68 

9.  80 

116.  75 

6  65 

25  36 

110.  75 

9. 

40 

27  00 

114 

62 

10.  75 

27.  04 

1C.  85 

117.  17 

7.  62 

25.  83 

111  58 

10 

46 

27  57 

115 

00 

11.  75 

27.  35 

11.  97 

117.  57 

8  55 

26  28 

112.  30 

11 

83 

28  06 

115 

38 

12  69 

27  69 

13.  10 

117.  79 

9.  47 

26  53 

112  70 

13 

39 

28  63 

115 

87 

13.  94 

28  00 

13.  90 

118.  00 

10  66 

2.  10 

113  30 

14 

85 

28  94 

116 

21 

15  08 

28.  28 

15  00 

118  13 

11  77 

27.  27 

113  79 

16 

46 

29  4  1 

116. 

42 

16.  08 

28.  55 

16  00 

118  25 

12.  81 

27  73 

114.  25 

19 

00 

29  76 

116 

70 

17.  12 

28  86 

17  05 

118.  42 

13  83 

28.  15 

114.  47 

21 

55 

29  98 

117 

79 

18.  19 

29.  09 

18  14 

118  75 

14  85 

28  41 

114.  79 

19.  37 

29.  36 

19.  30 

119.  04 

16.  45 

28  83 

115  14 

20.  43 

29.  63 

20.  46 

119.  21 

18.28 

29.  26 

115  66 

21.  61 

29.  85 

21.  75 

119.  34 

20.  15 

29.  50 

116.  79 

22.  70 

30.  07 

22.  74 

119.  50 

21.  50 

29.  63 

117.  58 

23.  50 

30.  20 

1 

Run  No,  7 
«  •  78  ft 
U.  -  20.25  fp» 

-  118.0  °F 

T  •  40. 0  °F 
0 

z  V  T  z 

in.  fpg _ °F _  in, 

0  125  8  85  83  87  0  125 

0.350  9  76  87.43  0  25 

0.50  11  50  90.87  0.50 

0.75  12  18  93.43  0.75 

1.0  12.68  95.22  1.0 

1.5  13.30  97.48  1.5 

2.0  14.00  100.04  2.0 

2.5  14.48  102.17  3.0 

3  0  14.91  103.34  4.0 

3.5  15.13  104.46  5.0 

4.0  15.38  105.58  6.0 

5.0  15.78  108.04  7.0 

6.0  16.09  109.21  8.0 

7.0  16.38  110.54  9.0 

8  0  16.79  111.38  10.0 

9.0  17.08  112.38  11.  C 

10.0  17.35  113.17  12.0 

11.0  17.84  113.70  i3.0 

12.0  17.81  114.34  14.0 

13.0  18.17  114.70  15.0 

14.0  18.38  114.96  16.0 

15.0  18.53  115.21  18.0 

16.0  18.79  115.51  20.0 

17.  0  18.  95  115.  79  22.0 

18.0  19.23  116.08  23.0 

19.0  19.47  116.38 

20.0  19.71  116.92 

21.0  19.87  117.25 

22.0  20.03  117.38 

23.0  20.07  117.79 

23.  75  20.  12  118.  87 


Run  No.  8 

Run  No.  9 

l  •  78 

ft 

X  • 

78  ft 

U„  .  19. 

95  fps 

u«  * 

20.15  fp« 

T„  »  116 

.2  OF 

T  ■ 

m 

111.0  °F 

T  •  40. 
0 

0  °F 

To  ■ 

37.0  °F 

U 

T 

Z 

U 

T 

fps 

Op 

in. 

fp« 

Op 

8.  77 

79.  92 

0  094 

8  58 

75  75 

10.  07 

83.  87 

0  16 

9  90 

80  70 

11  36 
11.  78 

87  83 

89  92 

■ 

10  66 
11.  64 

82  75 

86  75 

12  16 

92.  26 

1.  14 

12  28 

89  50 

12.  93 

95.  00 

1.  *Z 

12.  93 

92.  00 

13.  63 

96  78 

2.  13 

13.  53 

93.  00 

14.44 

100.  13 

2.  63 

13.  87 

95.  4 

!5.  00 

102.  96 

3  24 

14.  10 

96.  5 

15.  43 

105.  08 

3.  86 

14.  35 

98.  5 

15.  86 

106.  74 

4.  94 

14.91 

100  6 

16.  21 

108.  04 

6.  03 

15.  46 

103.  0 

16.  49 

109.  22 

7.  06 

15.  87 

104.  0 

18.  87 

119  35 

B.  10 

16.  17 

105.0 

18.  87 

111.  40 

9.  27 

16.  48 

106.  0 

17.15 

112.  35 

10.  37 

16.  85 

17.43 

113.04 

11.  50 

17.  23 

107.0 

17.  72 

113.  50 

12.  80 

17.60 

108.0 

18.  00 

114.  00 

14.20 

17.  97 

108.85 

18.  27 

114.  43 

15.  50 

18.39 

109.  0 

18.  44 

114.57 

16.  55 

18.73 

109.  5 

18.  86 

114.78 

17.80 

19.  16 

109.  8 

19.  28 

115.  26 

18.75 

19.38 

110.  04 

19.70 

115.  60 

20.  00 

19.  64 

no.  13 

19.  80 

116.  00 

21.06 

19.  82 

110.  22 

22.25 

19.  93 

110.  44 

23.  25 

20.  00 

110.  65 

Kim 

x  * 

7«S  t  ' 

11^  » 

JO. -0  Ips 

To  » 

117.5  °l 

ro  * 

40.0  °K 

z 

U 

in. 

fp» 

il 

0.125 

8.82 

0. 

0.  140 

9. 17 

0., 

0.34 

9.98 

0.. 

0.50 

10.44 

o.: 

0.74 

11.17 

0. 

0.96 

11.72 

i.i 

1.26 

12.24 

l. 

1.62 

12.85 

l. 

2  13 

13.47 

2. 

2.83 

13.95 

3. 

3.49 

14.34 

4. 

4.17 

14.75 

5. 

5.03 

IS.  13 

6. 

5.90 

15.52 

7. 

6.77 

15.83 

8. 

7.50 

16. 13 

9. 

8.34 

16.34 

10.. 

9.23 

16.59 

11. 

10.25 

16.96 

12. 

11.25 

17.32 

13. 

12.2S 

17.67 

14. 

13.50 

18.05 

15. 

14.70 

18.31 

16. 

1S.SS 

18.57 

18. 

lb. 60 

18.87 

19. 

18.01 

19.17 

20. 

18  65 

19.41 

21. 

19.73 

19.68 

22. 

20.73 

19.82 

21.77 

19.98 

22.75 

20.08 

23.50 

20.  IS 

'll  \N  VI I  *  H  I IV  \NI>  ll'll'l  IIUIIIII  ItATA  -  Continued 


■  111,  HO  fpn 
-  118.4 


108. 

92 

110, 

00 

111. 

08 

112. 

30 

113. 

00 

Run  No.  12 

*  •  70  ft 
•  19.80  fps 

-  118.4  °F 
T„  -  39.0  °F 


0.504 

11.46 

92.34 

0.804 

12.24 

95.61 

1.14 

12.84 

97.75 

1.51 

13.28 

99.70 

1. 88 

13.64 

101.00 

2.49 

14. 17 

103.57 

2.90 

14.40 

104.00 

3.48 

14.94 

103.8.3 

3.94 

15. 13 

106.79 

4.85 

15.67 

108.70 

S.  63 

15.99 

109. S3 

6.52 

16.35 

111.00 

7.62 

16.76 

112. SO 

8.65 

17.11 

113.00 

9.65 

17.49 

113.38 

10.79 

17.83 

113.79 

12.09 

18.  10 

114.21 

13.52 

18.51 

114.53 

15.22 

18.86 

115.04 

17. 3S 

19.29 

115.34 

19.23 

19.44 

115.86 

21. S6 

19.52 

117.38 

94 


TABLE  II 


MEAN  VELOCITY  AND  TEMPERATURE  DATA  -  Continued 


z 

in. 

Run  No .  15 

x  •  78  ft 

U_  -  9.70  fpj 

T  •  117.0  °F 

m 

T  ■  58.0  °F 

0 

U 

fp» 

T 

°F 

Run 

x  - 

u.  * 

T.  ■ 

To* 

Z 

in. 

No.  14 

78  ft 

9. 55  fps 

117.0  °F 

58.0  JP 

U 

fp* 

l 

in. 

Run  No.  15 

x  -  78  ft 

U„  •  10.10  fp» 

T.  •  117.0  °F 

T  •  58,0  °F 
o 

U  T 

fpi  °F 

0.062 

2.55 

72.6 

.062 

2.48 

0.125 

5.17 

75.78 

0.125 

5.48 

76.4 

0.125 

5.41 

0.25 

4.25 

79.96 

0.250 

5.88 

79.7 

0.250 

5.91 

O.SO 

4.76 

82.74 

0.  575 

4.16 

81.5 

0.575 

4.20 

0.75 

5.10 

85.15 

0.50 

4.57 

8'. 7 

0.50 

4.44 

1.0 

5.40 

87.09 

0.75 

4.66 

8S.6 

0.75 

4.03 

l.S 

5.82 

90.22 

1.0 

4.96 

87.6 

1.0 

5.10 

2.0 

6.17 

92.87 

1.5 

S.54 

90.8 

1.5 

5. SI 

2.5 

6.47 

95.00 

2.0 

5.75 

95.2 

2.0 

5.87 

5.0 

6.75 

97.04 

2.5 

5.98 

96.0 

x.  5 

6.09 

4.0 

7.07 

100. S2 

5.0 

6.29 

98.0 

5.0 

6.79 

S.O 

7.50 

102.91 

5.5 

6.52 

99.9 

5.5 

6.41 

6.0 

7.49 

104.96 

4.0 

6.69 

101.4 

4.0 

6.64 

7.0 

7.70 

106.47 

4.S 

6.91 

102.9 

4.5 

6.80 

8.0 

8.00 

107.58 

5.0 

7.11 

105.8 

5.0 

6.94 

9.0 

8.20 

108.58 

5. 5 

7.20 

104.7 

S.S 

7.05 

10.0 

8.54 

109.54 

6.0 

7.56 

105.7 

6.0 

7.17 

11.0 

8.50 

110.00 

7.0 

7.58 

107.2 

7.0 

7.41 

12  0 

8.68 

110.50 

B.O 

7.76 

108.0 

8.0 

7.65 

15.0 

8.77 

110.79 

9.0 

7.97 

109.1 

•*.0 

7.82 

14.0 

8.90 

111.21 

10.0 

8.14 

110.0 

10.0 

7.97 

15.0 

9.04 

l"  45 

11.0 

8.59 

110.4 

11.0 

8.15 

16.0 

9.15 

111.70 

12.0 

8.59 

110.9 

12.0 

8.27 

17.0 

9.25 

112.08 

15.0 

8.75 

111.5 

15.0 

8. 59 

18.0 

9.59 

112.47 

14.0 

8.86 

111.7 

14.0 

8.  SO 

19.0 

9.55 

112.92 

16.0 

9.09 

112.1 

15.0 

8.60 

20.0 

9.72 

115.50 

18.0 

9.52 

112.8 

l'-.0 

8.71 

21.0 

9.89 

115.94 

20.0 

9.49 

115.8 

17.0 

8.80 

22.0 

9.94 

114.83 

22.0 

9.S6 

115.9 

18.0 

8.89 

25.0 

9.99 

116  04 

19.0 

8  97 

25.75 

10.09 

117.04 

20.0 

9.04 

22.0 

9.18 

95 


in. 


0. 125 
0.  ISO 
0.46 
0.72 
0.99 
1.26 
1.61 
1.99 
2.65 
3.09 
3.S8 
4.02 
4.50 

4.98 
5.38 
5.89 
6.88 
8.03 
0.25 

9.98 
11.07 
12.00 

13.11 

14. 10 
’.4.96 
16.02 

17.12 

18.10 
19.  OS 
20.00 
21.15 
22.40 
23.50 


TABLE  II 

MEAN  VELOCITY  AND  TEMPERATURE  DATA  -  Continued 


No.  16 

Run  No .  17 

78  ft 

x  .  78  ft 

10.05  fps 

U.  *  in. 7  fpj 

114.0  °F 

T«  •  115.5  °F 

40.0  °F 

T  •  40.0  °F 

0 

U 

T 

z 

V 

T 

fpj 

Oc 

in. 

fp* 

°F 

3.21 

72.7 

0.084 

KS3 

4.21 

77.5 

0.114 

3.80 

4.74 

82.0 

0.  144 

4.10 

KX3 

5.13 

84.7 

0.254 

4.46 

81.30 

5.40 

36.3 

0.480 

4.97 

84.75 

-- 

88.7 

0.79 

5.26 

86.34 

5.78 

1.29 

5.70 

89.79 

6.  10 

92.4 

1.76 

6. 15 

92.50 

6.42 

94.7 

2.22 

6  35 

95.00 

6.64 

96.3 

2.79 

6.79 

97.66 

6.83 

97.6 

3.46 

7. 12 

100. 10 

7.00 

99.4 

4.21 

7.43 

102.30 

7. 13 

100.6 

S.  17 

7.76 

104.00 

7.25 

101.5 

6. 19 

8.05 

105.92 

7.  35 

102.0 

7.45 

8.42 

108.58 

7.58 

103.4 

9.03 

8.70 

109.43 

7.99 

104.5 

11.83 

9.05 

111.12 

8. 19 

■T7Y^H§Pf 

21.44 

10.62 

115. 'C 

8.36 

lot. 7 

8.47 

107. 1 

8.61 

107.9 

8.74 

-- 

8.88 

109.1 

8.97 

109.  S 

9. 11 

109.8 

9.23 

-- 

9.37 

9.46 

-- 

9.54 

-- 

9.62 

9.79 

111.0 

9.90 

112.3 

9.96 

113.7 

96 


TABuE  ill 


DATA  ON  RMS  TURBULENT  FLUCTUATIONS 


Run  No.  1 

Run  No.  4 

U  « 

no 

30.05  fps 

U^  ■  30.4  fps 

r  * 

m 

119.9  °F 

T  -  119.6  °F 

00 

T  « 

0 

39.0  °F 

T0  -  40.0  °F 

z 

v'FT  z 

z 

fix 

z 

in. 

Op 

in. 

fpS 

in. 

fps 

in. 

0.062 

5.48 

0.12S 

2.510 

0.13 

1.863 

SUB 

0.125 

S.  16 

0.34 

2.532 

0.21 

1.896 

mSm 

0.25 

4.89 

0.54 

2.464 

0,30 

1.908 

KB 

0.50 

4.62 

0,96 

2.361 

0.47 

1.898 

sum 

0.75 

4.44 

1.41 

2.320 

0.67 

1.812 

0.68 

1.0 

4.26 

2.0 

2.187 

1.03 

1.804 

1.01 

1.5 

3.0 

2.121 

1.43 

1.764 

1.39 

3.80 

2.07 

1.727 

1.87 

2.5 

3.58 

5.8 

1.883 

2.81 

1.705 

2.55 

3.35 

1.793 

3.56 

1.691 

3.45 

3.5 

3.17 

9.9 

1.655 

4.35 

1.687 

4.42 

4.0 

3.08 

11.9 

1.539 

S.04 

1.689 

5.00 

4.5 

2.94 

14.75 

1.470 

5.90 

1.689 

5.80 

EWiKIB 

2.76 

17.8 

1.350 

6.60 

1.665 

6.60 

5.5 

2.44 

7.60 

1.614 

7.70 

2.17 

8.65 

1.561 

8.70 

|sS|| 

1.99 

9.70 

1.552 

9.80 

1.77 

10.75 

1.538 

10.85 

1.54 

11.75 

1.513 

11.97 

1.27 

12.69 

1.494 

13.10 

13.94 

1.467 

13.90 

12.0 

0.91 

15.08 

1.425 

15.00 

0  72 

16.08 

1.403 

16.00 

17.12 

1.349 

17.05 

18.19 

1.290 

18.14 

18.0 

0.41 

19.37 

1.230 

19.30 

0.40 

20.43 

1.178 

20.46 

22.0 

0.39 

21.61 

1.116 

21.75 

22.70 

1.047 

2c.  74 

\j  w'2 

fps 


1,139 

1.253 

1.358 

1.375 
1.393 
1.475 
1.425 
1.441 
1.441 
1.455 
1.455 
1.437 
1.424 
1.409 
1.403 
1.395 
1.384 

1 . 376 

1.376 
1.369 
1.352 
1.319 
1.288 
1.263 
1.231 
1.190 
1.154 
1.077 
1.035 
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TABLF  III  || 

|  • 

DATA  ON  RMS  TURBULENT  FLUCTUATIONS  -  Continued  |  j 

1  i 


Run 

No.  7 

Run  No. 

10 

U  * 

00 

20.2a  fps 

U^  =  20.2  fps 

T  = 

CO 

118.0  °F 

T  =  117. 

yy 

5  °F 

To  * 

40.0  °F 

T  *  40.0  °F 

0 

z 

\l  V~2 

z 

z 

\f  v'2 

z 

\/w'  2 

in. 

°F 

in. 

fps 

in. 

fps 

in. 

fps 

0.062 

5.14 

0.25 

1.532 

0.125 

0.984 

0.125 

0.591 

0.12S 

4.73 

0.5 

1.483 

0.20 

1.038 

0.140 

0.634 

0.25 

4.47 

l.C 

1.374 

0.37 

1.012 

0.34 

0.753 

O.SO 

4.06 

1.5 

1.310 

0.57 

0.996 

0.50 

0 . 768 

0.75 

3.90 

2.0 

1.254 

0.74 

0.990 

0.74 

0.789 

1.0 

3.73 

3.0 

1.210 

1.02 

0.970 

0.96 

0.792 

1.5 

3.65 

4.0 

1.188 

1.39 

0.964 

1.26 

0.797 

2.0 

3.48 

5.0 

1.133 

1.90 

0.982 

1.62 

0.813 

2.5 

3.32 

6.0 

1.047 

2.58 

0.972 

2.13 

0.815 

3.0 

3.23 

7.0 

0.988 

3.26 

0.961 

2.83 

0.814 

3.5 

3.11 

9.0 

0.942 

3.85 

0.940 

3.49 

0.814 

4.0 

2.98 

11.0 

0.897 

4.39 

0.935 

4.17 

0.813 

4.5 

2.86 

13.0 

0.842 

5.20 

0.926 

5.03 

0.817 

5.0 

2.74 

15.0 

0.788 

6.15 

0. 970 

5.90 

0.803 

5.5 

2.37 

18.0 

0.713 

7.20 

0.922 

6.77 

0.807 

6.0 

2.40 

21.0 

0.554 

8.18 

0.921 

7.50 

0.813 

7.0 

2.16 

9.33 

0.922 

8.34 

0.809 

8.0 

1.87 

10.48 

0.922 

9.23 

0.817 

9.0 

1.70 

11.26 

0.923 

10.25 

0.820 

10.0 

1.49 

12.37 

0.893 

11.25 

0.811 

11.0 

1.33 

13.67 

0.8711 

12.25 

0.807 

12.0 

1.16 

14.73 

0.850 

13.50 

0.808 

13.0 

0.99 

15.79 

0.815 

14.70 

0.795 

14.0 

0.87 

16.80 

0.780 

15.55 

0.770 

16.0 

0.61 

18. 12 

0.760 

16.60 

0.748 

18.0 

0.46 

19.13 

0.703 

18.01 

0.720 

20.0 

0.44 

20.34 

0.658 

18.65 

0.700 

22.0 

0.41 

21.49 

0.622 

19.73 

0.668 

22.73 

0.578 

20.73 

0.646 

23.50 

0.500 

21.77 

0.616 

22,75 

0.579 

_ 

23.50 

0.520 

4 
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TABLE  III  1 


DATA  ON  RMS  TURBULENT  FLUCTUATIONS  -  Continued 


N  B 

_ ll 

Run 

U  «  V 

T  =  1 
<** 

T  -  3 
o 

to.  15 

0.10  fps 

17.0  °F 

8.0  °F 

2 

in. 

Run  No 

U  =  10. 

m 

T  =  114 

e» 

T  =  40. 
o 

.  16 

OS  fps 

.0  °F 

0  °F 

z 

in. 

i 

« 

i 

V  ^2 

Or 

\l  u'2  '  2 

*  1 

fps  |  in. 

\i~ 

fps 

fps 

0.12.- 

4.72 

0.125 

1.772 

0.125 

0.148 

0.125 

0.124 

1  0.2S 

4. 4c 

0.23 

0.641 

0.’5 

0.259 

0.19 

0.20?  | 

0.50 

3.94 

0.50 

0.561 

0.46 

0,299 

0.55 

0. 262  * 

i.O 

3.55 

1.0 

0.553 

0.72 

0.304 

0.88 

0.27?  ! 

1.5 

3.52 

1.5 

0.526 

0.99 

0.312 

1.22 

0.27b  i 

2.0 

5.44 

2.0 

0.511 

1.26 

0.316 

1.71 

0.274  i 

2.S 

3.52 

3.0 

0.495 

1 . 61 

0.318 

2.37 

0.274  j 

3.0 

3.15 

4.0 

0.47.3 

1.99 

0.318 

2.85 

0.269  1 

3.5 

2.98 

5.0 

0.407 

2.65 

0.325 

3.31 

0.266  | 

4.0 

2.74 

6.0 

0.392 

3.09 

0.329 

4.26 

0.261  * 

4.5 

2.49 

7.0 

0.382 

7.02 

0.336 

5.50 

0.263  ! 

5.0 

2.20 

9.0 

0.336 

4.50 

0.332 

6.60 

0.268  ! 

6.0 

1.87 

31.0 

0.308 

4.98 

0.528 

8.3C 

0.273  J 

7.0 

1.70 

13.0 

0.289 

5.38 

0.323 

9.73 

0.272 

8.0 

1.41 

15.0 

0.289 

5.89 

0.321 

11.15 

0.275  S 

9.0 

1.20 

18.0 

0.290 

6.88 

0.324 

12.72 

0.276 

10.0 

0.99 

8.03 

0.332 

14.40 

0.277 

11.0 

0.95 

9.25 

0.331 

16.35 

0.Z72 

12.0 

0.61 

9.98 

0.330 

18.20 

0.245 

13.0 

0.50 

12.00 

0.323 

14.0 

0.44 

14.10 

0.313 

> 

16.0 

0.40 

16.02 

0.297 

18.0 

0.35 

13. 10 

0.253 

20.0 

0.31 

22.0 

0.27 

» 

t 
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DATA  ON  ONE -DIMENSIONAL  SPECTRA  -  Continued 
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TABLF  V 


DATA  ON  ONE -DIMENSIONAL  SPECTRA  -  Continued 
U_  »  20  fpx  ,  R  -  0.0254 

*J| 


i  •  0. 

ki  ♦» 

ft-1 

125  in. 
.(kjJ/t7? 
ft 

t  -  O.SO 

*1  *t 

ft'1 

in. 

(k,j/r-T 

ft 

t  -  2.0 

kl  *t 

ft'1 

in. 

(kiJ/t^ 

ft 

:  ■  5. 

ki 

ft*1 

0  in. 

ft 

■ 

«N 

m 

1.13x10* 

1.26x10- 2 

9.14 
l.  14x10* 

1.60xl0"2 

7.50 

2. 24xlO*2 

6.72 

2. SlxlO*2 

6.00 

3.caio-J 

1.41 

1.34 

1.66 

9.38 

2.55 

8.40 

3.00 

7.51 

3.19 

1.76 

1.23 

1.43 

1.67 

1.17x10* 

2.16 

l.OSxlO1 

2.75 

9.39 

2.81 

2.26 

1.19 

1.83 

1.48 

1.50 

2.18 

1.34 

2.60 

1.20X101 

2.81 

2.82 

1.02 

2.29 

1.39 

1.88 

1.78 

1.68 

2.14 

1.50 

2.40 

3.55 

9.3Sxl0'3 

2.86 

1.01 

2.34 

1.32 

2.10 

1.64 

1.88 

1.66 

4.S2 

7.S6 

3.66 

8.S2X10-3 

3.00 

1.03 

2.69 

1.36 

2.41 

1.21 

5. 65 

5.72 

4.57 

6.S4 

3.75 

S.OSxlO"3 

3.36 

9.47X10*3 

3.00 

8.46x10*® 

7.06 

3.89 

5.71 

4.71 

4.69 

5.88 

4.20 

6.47 

3.76 

6.12 

8.82 

2.90 

7.15 

5.S7 

5.86 

4.39 

5.25 

4.46 

4.69 

4.68 

l.lixlO2 

2.25 

915  ? 

2.62 

7.  SO 

2.97 

6.72 

3.20 

6.00 

3.09 

1.41 

i.24 

1.14XHT 

1.80 

9.38 

i.79 

8.40 

2.14 

7.51 

1.81 

1.76 

S. 99x10 

1.43 

9.18X10'4 

1.17xl02 

9.2OX10-4 

l.OSxlO2 

1  - 1 2 

9.39 

9.63xl0'4 

2.26 

3.91 

1.83 

6.23 

1.50 

5.96 

1.34 

7.07xl0’4 

1.20x10* 

7.31 

2.82 

2.09 

2.29 

3.«7 

1.88 

3.44 

1.68 

3.40 

l.SO 

3.72 

3. S3 

1.15 

2.86 

1.05 

2.34 

1.91 

2.10 

2.03 

1.88 

1.99 

4. 52 

3.65xl0-S 

3.66 

8.32xl0-s 

3.00 

7.37x10-S 

2.69 

8.09x10'® 

2.41 

6.65X10"3 

5.65 

1.41 

4.S7 

3.16 

3.75 

3  23  X 

3.36 

2.95 

3.00 

2.89 

7.06 

4.67xi06 

5.71 

IIS  . 

4.69 

8.5SX10*6 

4.20 

9.95x10-* 

3.76 

1.00 

8.82 

1.22 

7.15 

3.36x10"® 

5.86 

1.32 

5.25 

3.44 

4.69 

3.37xlO'7 

1.13xl03 

4. 87x10" 7 

9.  IS 

l.OS 

7.50 

6.22x10-' 

6.72 

7.40xl0‘7 

6.00 

7. 39x10" 7 

1.41 

1.76 

2.76 

1.77 

1. 14xl03 
1.43 

3.59xlC* ' 

i.:s 

9.38 

1.42 

6.40 

l.OSxlO 

1.73 

5.66x10'* 

7.51 

2.14 

U_  •  10  fps  .  R,  -  0.113 

«  i 


t  -  0 

12S 

X  •  O.SO 

in. 

x  -  2.0 

in. 

2-5 

0  in. 

X  -  10. 

0  in. 

kl  ♦, 

kl  *V 

(kjJ/t7* 

kl  4t 

(k^/t7! 

kl 

♦t,(ki)/t^ 

kl 

^.(ki)/!7? 

ft-1 

ft 

ft*1 

ft 

ft’1 

ft 

ft-1 

ft 

ft'1 

ft 

3.09X101 

1.2SX10'2 

2.09x10* 

1.28xl0'2 

1.66X101 

1./4X10'2 

1. 37x10* 

1.99xl0'2 

1.19x10* 

3.46X10-2 

3.87 

1.24 

2.62 

1.31 

2.07 

1.75 

1.72 

2.14 

1.48 

2.70 

4.83 

1.02 

3.27 

1.28 

2.S9 

1.53 

2.14 

1.91 

1.85 

1.93 

6.18 

8.91x10"® 

4.19 

123 

3.32 

1.53 

2.75 

1.91 

2.37 

1.92 

7.73 

5.35 

5.23 

8.02X10-3 

4. IS 

116  T 

3.43 

1.35 

2.96 

1.  IS 

9.67 

3.93 

6.54 

6.62 

5.19 

9.04x10*® 

4.29 

1.02 

3.70 

8.30X10-3 

1.24xl02 

197  4 

8.37 

3.77 

6.64 

S.4S 

5.49 

6.94xl0-3 

4.74 

6.83 

1.55 

8.68x10' 

1. 05x10' 

2.17 

8.29 

3.10 

6.86 

3.83 

S.93 

3.10 

1.93 

3.94 

1.31 

102  . 

1.04x10* 

1.66 

8.58 

1.98 

7.40 

1.81 

2.42 

1.73 

1.64 

S.13X10'4 

1.30 

8. 73x10“  4 

1.07xl02 

8.63xl0-4 

9.26 

9.97xl0"4 

3.09 

5.96x10"'* 

2.09 

2.25 

1.66 

3.78 

1.37 

4.58 

1.19x10' 

4.67 

3.87 

9.53x10'® 

2.62 

7. 23xl0-s 

2.07 

1.28 

1.72 

1.78 

1.48 

1.84 

4.83 

2.90 

3.77 

229  6 

2.59 

2.98X10’5 

2.14 

4.70x10° 

1.85 

5.  15x10"5 

6.18 

i.n 

4.19 

5.47x10** 

3.32 

1.05 

2.75 

1.51 

2.37 

1.46 

7.73 

3.22x10 
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Fig.  2  Probes  mounted  on  the  carriage 


Fig.  3  Instruments:  1.  Spectrum  analyzer  2.  Oscilloscope 
3.  Transonic  pressure  meter  4.  Hot-wire  anemometer 
5.  rms  volt-meter 


Figure  4.  Resistance  thermometer:  circuit  diagram 


Definition  sketch  of  momentum  and  thermal 
boundary  layers 


Figure  10.  Comparison  of  the  velocity  profile  at  x  *  78  ft 
with  that  at  x  *  70  ft,  *  20  fps 


z,  inches 

Figure  11.  Coaparison  of  the  temperature  profile  at  x  =  78  ft 
with  that  at  x  =  70  ft,  Ux  =  20  fps 
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Figure  12.  Variation  of  R,  with  z/5  ;  U  ■  30.0  ,  19.8 
and  10.7  fps 


distributions 


Figure  14,  Mean  temperature  distributions  in  the  boundary 
layer;  R,  «  0.0111  ,  0,0262  and  0,113 
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Figure  15.  Distribution  of  turbulent  intensities; 

R  =  3.02  x  105  and  R.  =  0.0112 
e5  *6 


:-C*ppe«* 


Figure  17.  Distribution  of  turbulent  intensities; 


Distribution  of  turbulent  energy  in  the  boundary 
layer;  R.  =  0,0112  ,  0.02S  and  0.108 


Distribution  of  viscous,  turbulent  and  total 
shear  stresses;  U  =  10.7  fps  ,  R.  *  0.0883 


*0.0883 
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Figure  24.  Distribution  of  horizontal  heat  flux; 
R  •  0.0113  ,  0.0254  and  0.104 
*6 


Distribution  of  the  ratio  -u't'/V't 
R.  =  0.0113  ,  0.0234  and  0.104 
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Distribution  of  momentum  exchange  coefficient; 
R.  =  0.0112  ,  0.0247  and  0.0883 
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Figure  30.  Distribution  of  the  ratio  of  exchange  coefficients  of  heat  and 
momentum;  R.  =  0.0112  ,  0.0247  and  0.0883 


Figure  31.  Comparison  of  mean  velocity  distribution  in  the 

wall  layer  with  Monin  and  Obukhov’s  similarity  theory 


with  similarity  theory 


Figure  33.  Variation  of  dimensionless  shear  coefficient 
with  stability  ratio 


30.05 

3000 

20.25 


Figure  37a.  Vacation  of  K^/K^  as  determined  from 
measured  fluxes  with  z/L 


Figure  37b.  K^/K^  as  determined  from  the  ratio  S/R 


function  of  z/L 


atmospheric  data 
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Figure  42.  y  w’  /u*  vs  Ri  :  comparison  with 
atmospheric  data 


z/S  =001 


comparison  with 


Figure  44.  Variation  of  y  with 

w  t 
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figure  48.  -<opres  ntation  of  w’-  spectra  in  the 
similarity  coordinates  and  comparison 
with  Heisenberg's  theoretical  expression 
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Figure  49.  Normalized  spectra  of  temperature 
fluctuations;  =  20  fps 
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